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The thesis presents a study of proton spin-lattice 
relaxation in gas mixtures H 2 “Ar, II 2 ~N 2 and H 2 -N 2 0, T^ 
minimum study in pure CH^Cl and an analysis of the existing 
data in NH^ P £ and HC1. 

Proton spin-lattice relaxation times T^ were measured 
as a function of density and temperature in hydrogen mixtures. 
This data together with the measurements done by Poster et al. 
in H 2 -Ar mixtures and by Williams in II 2 -N 2 mixtures were 
analysed using Bloom-Oppenheim theory to extract information 
on the anisotropic intermolecular interaction between these 
molecules. In CH^Cl, T-^ were measured as a function of 
density in the low density region at room temperature and 
T-j^ minimum was obtained. Prom the experimental value of T-^ 
minimum the effective spin-rotation coupling constant for 



protons in CliUCl was evaluated. The data obtained by- 
Courtney and Armstrong in By Tward and Armstrong in HC1 

and by Armstrong and Hanrahan in NH^ and BF ^ were analysed 
using Bloom-Oppenheim theory and a model first proposed by 
Rajan et al. for finding out the correlation time for spin- 
rotation interaction was tested. 

All the measurements were done with a spin-echo spec- 
trometer having an operating frequency of 30 MHz using phase 
coherent detection. The subject matter is divided into six 
chapters as described below; 

The Chapter I contains an introduction to the basic 
principles of Nuclear Magnetic Resonance and the phenomenon 
of spin-lattice relaxation. 

In Chapter II the method of spin-echoes to measure the 
relaxation times and the spin-echo spectrometer used for this 
purpose are ue scribed together with the details of the sample- 
holder and sample preparation. 

The experimental data and its analysis in the mixtures 
of H 2 with Ar, ]$2 ^0 a rev ^- ew the Bloom-Oppenheim 

theory are presented in Chapter III. T^ was measured as a 
function of density (10 < p<70 amagats) and composition in 
the temperature range of 300-600K for H^-Ar and H 2 -N 2 and 
200-400K for mixtures. Bloom-Oppenheim two level 

theory for proton spin-lattice relaxation in H 2 predicts a 



non-linear behaviour of T-^/p as a function of composition of 
the mixture PL,-X where X may be Ar, or 1^0 . • Since the 
fractional population of the rotational states higher than 
J=3 for ortho- end J=4 for para-Iig are not negligible at 
temperatures above 45 OK, this theor;*- can not be applied at 
higher temperatures. In the absence of an exact theory which 
is applicable to higher temperatures, the T-^/P was 

fitted to the following empirical relation using the method 
of least squares as was done by some previous workers for 
Hg-He mixtures 

Vp = Ao + + vl 


where A Q , A^ and A2 are constants and is the fractional 

concentration of X in the mixture. (T^/p) for 2^=1 gives 

the value of i^/p due to H 2 ~X collisions alone and is 

denoted by ( T^/ p ) Since the deviation from linear 

dependence is not very prominent up to 63$ of the solvent 

gas, the data was also fitted with a linear fit to obtain 

the value of ( f^/p) The ( T-^/P ) ^ A values obtained by 

the empirical fit and a linear fit differ by 8$ at 300K and 

by 2$ around 600K. Model potentials for the intermolecular 

interactions were proposed in all the three cases and the 

strengths of the attractive and repulsive terms were obtained 

0 **“ 

using the experimental value of (l\/P) . It was found 

that the values of the relative anisotropy in the attractive 



xii 

term of the intermolecular potentials for Ilg-Ar and 
are 0.137 and 0.131 respectively whereas the polarizability 
data gives a value of 0.128. Agreement between the values 
obtained from this analysis ana from polarizability data 
justifies the use of the proposed model potential to explain 
the UMR data in these systems, for the quadrupole 

moment of calculated from the strength of the attractive 

2 c 2 

term was found to be 4.65 x 10" esu cm while the recommended 

— 26 2 

value is 3.0 x 10 esu cm . The potential proposed for this 
system may be taken as a tentative potential. 

In Chapter IV the measurements in CH^Cl have been 
reported. T^ were measured in Crl^Cl gas in low density region 
(0.2 < p<6.5 amagats) at room temperature. Since T^/p is 
very long (370 msec/amagat) a pulse train consisting of tc/2- 
burst followed by y:/2-~n pulses was used to save the experimen- 
tation time. the spectrometer was modified for this purpose. 
From the T^ minimum the value of the effective spin-rotation 
coupling constant, C eff , for proton in CEUC1 was obtained 
which is 0.32 KHz. 

Chapter V contains an application of the model first 
proposed and used by Raj an et al. to evaluate the correlation 
time of spin-rotation interaction in some gases. In this 
model, it is assumed that the correlation time of the spin- 
rotation interaction can be approximated by the life time of 
the molecule in a J state if spin-rot arion interaction is the 
dominant mechanism of relaxation in the gas. Using this 
model results (available in literature) for NH^, BF^, F ? 
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and HC1 have been analysed using Bloom-0 ppenheim theory to 
find the dominant multipole moments of these molecules. It 
was found that the quadrupole moment Q of Pg . an ^ dipole 
moment u of HC1 are comparable to the results obtained from 
other experiments while Q for BP~ and p for are quite 
different from the reported values. Thus it is concluded 
that this moded. works well in some cases, especially for sph- 
erical top molecules such as CH^ as discussed by Rajan et al. 
and also for linear molecules as is found in the present work. 
Perhaps the model does not work well in symmetric top mole- 
cules. Keeping in view the remarks made by Armstrong and 
Courtney about the possible inadequacy of the existing 
theory for spin-relaxation in symmetric top molecules it can 
be said that the discrepancy observed in the values of 
multipole moments in the present analysis may be either due 
to the inadequacy of the model proposed by Itajan et al or 
due to the inadequacy of the spin-relaxation theory for 
symmetric top molecules itself. 

Conclusions of the work presented are given an Chapter 
VI. The circuit diagrams of the different parts of the 
spectrometer are given in Appendix. 


# 



CHAPTER I 


INTRODUCTION 


I . 1 Outline of the Principle of NHR 

Magnetic resonance is a phenomenon found in systems 
that possess magnetic moment as well as angular momentum V Y 
and occurs in the presence of external static and oscillat- 
ing magnetic fields applied perpendicular to each other. 

The principles of magnetic resonance are given in detail 
1 2 

elsewhere * and are presented here in brief. 

A static magnetic field H Q applied along, say, Z 
direction produces an interaction Hamiltonian, which for a 
nucleus with spin angular momentum I and magnetic moment 
ft = where y is the gyromagnetic ratio, is given by 


It = -5-ir 0 = -yh 0 *i z 

(1.1) 

The allowed energy levels are 


E m = H o m 

(1.2) 

where m = -I , -1+1 , . . . , +1 . 



The fractional populations of these energy levels at equili- 
brium are given by 



2 


exp (-E /kT) 

p m = i 7 r- (1-3) 

m' ez P (-V kT) 

where k is the Boltzmann constant and T is the absolute 
temperature of the system, for 11 weakly interacting spins, 
the equilibrium magnetization is given by 


= H I- P m y ft m 
m 

(1.4) 

= N y 2 ft 2 l(I+l)/3kT 

(1.5) 


where the high temperature approximation yhH 0 /kT« 1 is 
used. 


The approach to equilibrium of a spin system disturbed 
from its thermal equilibrium is often describable by the 
phenomenological Bloch equations 


dl? + V + 

If = y fix H - X_ 


M 0 - M . 
+ -X k 


T 


1 


( 1 . 6 ) 


where T-^ and T£ are called the longitudinal and transverse 
relaxation times respectively and H is the magnetic field. 
If the relaxation effects are neglected, the equation of 
motion for ii in a reference frame rotating about Z axis 
with angular frequency w may be written as 

f§ = ytf x (H + &) (1.7) 

Prom this it is seen that is stationary in the rotating 
frame for H = H Q k If J - -yB^k, which implies that, in the 



3 


laboratory frame the magnetic moment ft precesses about the 

a . 

field H Q k -with an angular frequencj- & 0 = -yH q £. This 
frequency is called larmor frequency. 

The energy levels (eq.(l.2)) of a magnetic system, say 

/N 

nucleus, in an external static magnetic field H 0 k are 
discrete and can be detected by some form of spectral 
absorption if the system is subjected to an interaction 
that can cause transitions between them. In order to cause 
transitions between two levels having energy E^ and E 2 
(E-^<* E 2 ) » the interaction must have a nonvanishing matrix 
element joining them and should satisfy the conservation of 
energy, i.e. the interaction must be time dependent and of 
an angular frequency to such that 

ha) = E 2 - E-l (1.8) 

These requirements are very well satisfied by an 
alternating magnetic field of frequency w applied perpen- 
dicular to the static field. If the alternating field is 
written as 

ft-^t) = H-j- (i cos cat + j sin w t) (1.9) 

then $ = ft 1 (t) + H q £ (1.10) 

and a.. perturbation term is introduced in the total 


Hamiltonian 
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lit = -yh H-, (l_ cos wt + I sin w.t) (1.11) 

j. jc y 

9 ? ■ 

The matrix elements of ft "between the states m ana m' ,v = ; 
vanish unless m' = m + 1. Therefore, the transitions occur 
between levels adjacent in energy and 

* o> = E 2 - 3 1 = yfiH 0 (1.12) 

or to = y H q (1.13) 

Therefore, in order to produce transitions, the alternating 
field must have a frequency equal to the precession frequency 
of the magnetic system in the static field. 

With both the fields applied together, the equation 

-V 

of motion of M in the rotating frame is given by 


g-f = yt x [(H 0 +<£) £ + i] (1.14) 

= yM x H e ££ (1.15) 

where H eff = (H Q + ~) ’k + &,_ i (1.16) 

The eq. (1.15) tells that M precesses about H e ^^. with ah 


angular frequency -yHg^^. When s i and the 

magnetization, therefore, precesses about A magnetic 

moment that is parallel to the static field initially will 
then process in the YZ plane always remaining perpendicular 
to and, thus, periodically line up opposite to S Q . AH 
the energy it takes to tilt M away from H Q is returned in a 
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complete cycle. There is no net absorption of energy from 
the alternating field but rather alternately receiving and 
returning of energy. This phenomenon is known aS Nuclear 
Magnetic Resonance. 

At resonance the precession frequency is yH^ and 
therefore if R^ is switched on for a time , ft will precess 
about through an angle 0 = yB^t^. 0 can be made 180°(7t) 
or 90°(-n;/2) by choosing H-. and t appropriately. After ) " 

-1 vAJ J? | J 

^ i 

application of a tc/2 pulse M precesses in the XT plane t 

V. d 

about H q and induce s an ,emf in the sample coil kept perpen-^ ( 
dicular to H Q . This induced signal is known as Free 
Induction Decay (F.I.D.). If a % pulse is applied at a 
time t = t after the n/2 pulse, the spins tend to re phase 
in the XI plane at a time 2t resulting in a free induction 
signal known as a spin-echo. 


If the total magnetic field is written as 


$ = H 0 k + (2 Bj 


cos a) t) i , K, H 0 


(1.17) 


then the solution of the Bloch equations (eq. 1.6), in the 
absence of H^, are obtained as 


MjjXf) = M^O) [cos(a)t +4,)] exp (-t/T 2 ) (1.18) 

My( t) = \ 7 (0) [ sin(w t +<*>)] exp (-t/lk,) (1.19) 

M 2 (t) = M Q + [M z (0) - M q ] exp (-t/1^) 


( 1 . 20 ) 
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where +* My ^ and the relaxation effects have also 

been considered. 

By observing M^t) (or M (t)) and M z (t) as a f motion 
of time, the values of T „ and T^ respectively can be deter- 
mined experimentally. the measurements and interpretation 
of only T-^ are presented in this thesis. 

I . 2 A Note on Relaxation Mechanisms 

After the spin-system has been disturbed from equili- 
brium by an external agency, it tends to approach the thermal 
equilibrium again by transferring its excess energy to the 
other degrees of freedom (called 'lattice') . The process of 
this transfer of energy is called spin-lattice relaxation 
and occurs through different types of interactions within 
the molecule and outside. The interactions usually assumed 
to be responsible for thermal relaxation are: 

( a) The In t ramol e cular dipolar and quadrupolar (for I » 
only) interactions between neighbouring nuclei of the same 
molecule. These interactions transform as ft) under 

rotations of the molecule where ft denotes the orientation 
of a vector fixed to the molecule with respect to the spac 
fixed frame. 

(b) The Intermole cular dipolar interactions between nuclei on 
two different molecules. 
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(c) The in t ramo le cul ar spin - rotational interactions between 
nuclear spin and the rotational angular momentum of the 
molecules. The dominant terms of this interaction transform 
as Jj the rotational angular momentum of the molecule. 

It is usually assumed in the conventional theories 
for relaxation that all the three mechanisms mentioned above 
contribute independently to the relaxation process, so that 
one can write 

jjr- = R^ + Itg + (1.21) 

where R^, Rg and R^ are the contributions to the relaxation 

rate due to the mechanisms a, b and c respectively. In 

molecular gases which are not too dense, R^ is negligible, 

so that the relaxation rate ~~ can be written as 

i l 

T~ = il A + (1.22) 

Expressions for relaxation rates for H 2 were obtained 

3 4 

by Schwinger , Reedier and Opechowski and Johnson and 

5 

Waugh in terms of the correlation times of the intramolecular 
interactions assuming that only J=1 state is populated. 

Bloom and Oppenheim “ have developed a theory in which 
this restriction is removed and transitions between different 
J states are allowed. They have also obtained expressions 
for the correlation times of the intramolecular interactions 
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in terms of the intermole cular interaction by treating the 

dynamics of the system with ' Constant Acceleration Approxi- 

11 12 

Elation' (CAA). Kinsey et al. * Poster and Rugheimer 

13 

and Shafer and Gordon have done Quantum Scattering 

Calculations. Gordon^ - ^ has developed a theory based on 

kinetic theory of gases and obtained expressions for spin-- 

relaxation rates in terms of the intermole cular potentials 

for molecules in gas phase. For heavier molecules, where 

large number of rotational states are populated, Raj an 
17 

et al. have obtained expressions for 1'^ in terms of the 
intermole cular potentials assuming that the correlation 
time of the dominant intramolecular interactions is equal 
to the life time of the molecule in a J state. 

I. 3 Information from Tq Measurements 

( a) Intermole cular Potentials 

In a gas, the molecules are always colliding with 
each other in a random way. During the collisions the 
intermole cular forces cause reorientation of the molecules 
(i.e. change the rotational states) resulting in a change 
of the local fields set at the sites of the nuclei. Due 
to these randomly fluctuating local fields the nuclei 
approach thermal equilibrium. The rate, 1/i 1 -^, at which 
the nuclear spins approach the state of equilibrium can 
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be measured experimentally. Therefore, measurement of 
spin-lattice relaxation time T^ provides detailed information 
on the process of molecular reorientation which in turn 
gives information about the anisotropic interaction between 
the molecules. The distance of closest approach for the 
molecules decreases as the temperature increases and, thus, 
it is possible to probe more into the short range part of 
the intermole cular potential by measuring T^. at higher 
tempe ratures . 

There have been several previous experimental investi- 
gations of T-j, in Ii 2 10, 18-20 and in the mixtures of H 2 with 
other gases^ ,l2 ’ where information on the anisotropic 
part of the in ter mole cular potential has been obtained from 
the analysis of the proton T^ measurements of .the ortho-H 2 
molecules as a function of density, temperature and concen- 
tration of H 2 in the binary gas mixtures. However, not much 
work has been done at higher temperatures (1 > 300 K) . 

A part of the work presented in this thesis covers a 
study of proton spin-lattice relaxation in the mixtures of 
H. with mono-, di- and tri-atomic gases, namely, Ar, N ? 
and N 2 0 as a function of density, composition and tempera- 
ture. T^ were measured in H 2 ~Ar and H 2 -H" 2 mixtures in the 
temperature range of 300-600K whereas in H 2 -ET 2 0 in the 
range of 200-400K. These data together with ^ values 
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12 24 ' 

available for H^-Ar and mixtures below room 

temperature were analysed using Bloom-Oppenheim theory to 
extract information about the intermole cular anisotropic 
potentials. 


( b) Intramole cular Interactions 


The rate of spin-lattice relaxation is determined 
essentially by the component, at the 1 armor frequency, of 
the fluctuating fields at the site of nuclei ana depends 
upon the duration as well as on the strength of the inter- 
action responsible for relaxation. If the average time 
between collisions that change the state of the molecule 
be t , it is expected that ^ « t if the fluctuations are 

O C 

1 

very fast ( — » m)j i.e. shorter the % is, the smaller 
c 

is the component of the field at Larmor frequency. As t 
increases, the fluctuating fields will have maximum effect 


when 


1 


is of the order of thus goes through a 


minimum at 


1 


w Q . for a dilute gas it can be shown that 

-1 


at constant temperature r « p where p is the density 

c 

of the gas. Therefore, T^ as a function of density shows 
a minimum which is different for different temperatures. 

Since stronger the strength of the interaction between the 
nuclear spins and the fluctuating field, faster the 
rate of relaxation, the value of 1’^ minimum yields information 
on the strength of the intramolecular interactions. 
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One of the contributions of this thesis is the T^ 
minimum study in CH^Cl gas. 

I • 4 Scope of the Present Work 

The apparatus and experimental techniques that were 
used to obtain the T^ data are presented in Chapter II. In 
Chapter III the theory used to interpret the T-^ data in 
H^-Ar, H 2"’'^2 t ' ) mirfcures along with the analysis of 

the data to obtain information about the intermolecular 
potential is given. T^ minimum studies and the analysis of 
the data to obtain the value of the effective spin-rot a'c ion 
coupling constant for protons in CH^Cl are described in 
Chapter IY. Chapter V contains an application of a model 
first proposed and used by Rajan et al. 1 ^ to evaluate the 
correlation time of the dominant intramolecular interaction 
in molecules having large moment of inertia. Using their 
model T^ resin ts (available in literature) for ^ , 

2 ^ 07 

and HCx have been analysed to find the 
strengths of the dominant multi pole moments of these mole- 
cules. Chapter VI presents the conclusions. The circuit 
diagrams of various parts of the spectrometer are given 
in Appendix. 



CHAPTER II 


EXPERIMENTAL TECHNIQUES AHu APPARATUS 


I I . 1 The Method o f Spin -Echoes 

The nuclear spin relaxation times can be directly 

OQ 

measured using the method of spin-echoes ° which is described 
below; 

Let us consider a bulk sample containing large number 
of nuclear spin moments placed between the pole pieces of a 
strong magnet. When the field H Q , say in the Z~direction, 
is switched on, the tiny nuclear moments tend to align 
themselves along H Q giving rise to a state of thermodynamic 
equilibrium for the sample with a net macroscopic magnetiza- 
tion M q along the Z axis. As described in chapter I, if a 
radio frequency field H-, at Larmor frequency yH 0 , where y 
is the gyromagnetic ratio of the spins, is applied to the 
sample for a time , and in a direction perpendicular to 
H , the magnetization is tilted away from the Z direc- 
tion by an angle 0 given by 


6 


riy t u 


( 2 . 1 ) 
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Intense rf pulses are used to tip the magnetization M Q 

through an angle 0 in a time much shorter than or r £ so 

that the relaxation effects while is on can he neglected. 

In the absence of 13^, M q precesses about H Q . A m/2 rf 

pulse applied to the sample coil whose axis is along the X 

direction rotates the magnetization^ to the XY plane. let 

us visualize this situation in a frame X'Y'Z' rotating with 

rf and where Z' and Z coincide. Following the n/2 pulse 

applied along positive X' direction, the magnetization TVI^ 

2 o 

lies entirely along the Y' axis (Fig. la) J and sees the 

HK 

field H . Thus, when looked at from the laboratory frame 
XYZ, it precesses about H Q but lying in the XY plane . It 
induces an emf in the sample coil which can be detected and 
observed. This signal is known as 'Free Induction Decay' 
(FIi>). JJue to transverse relaxation this signal decays. 

The time constant of the decay in the X-Y plane would be T^ 
if the field H Q is perfectly homogeneous. However, for all 
the practical sample sizes, which cannot be too small due to 
signal amplification limitations, there exists an inhomo- 
geneity, say AH q , over the sample dimensions. So the nuclear 
moments, in the different parts of the sample experience 
different fields and hence precess with different Larmor 
frequencies. Thus, there is a spread of precession frequen- 
cies centered at oj q . If this situation is visualized from 
a reference frame rotating with frequency u 0 about z axis, 
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the spins will appear to fan out since some nuclei precess 
faster and some slower than the frame (Pig- lb). Thus the 
magnetization in X-Y plane decays in time of the order of 
(l/yAH 0 ) due to field inhomogeneity. Application of a n 
pulse at time t much smaller than or T^ after the n/2 
pulse rotates the fanning out moments through 180° about the 
X axis (Pig. lc) . So the phase differences produced 
following the n/2 pulse, due to field inhomogeneity, are 
reversed (Pig. Id). The spins continue to precess in the same 
direction as they were before the application of the a pulse. 
This leads to reclustering of the tiny spin moments (Pig. le) 
at a time t = 2% after the n/2 pulse and a maximum signal is 
built up in the coil which decays again due to same reasons 
given above. This second signal is known as ’spin-Echo’. 

Though the free induction decay can also be used to 
measure T^, the spin-echo has got, from experimental point 
of view, several advantages few of which are: (i) the n and 
n/2 pulses are off during the time the echo is observed. 

The observed signal is, therefore, free from the interference 
due to the pulses, (ii) when the echo appears the receiving 
system has already recovered from saturation and other effects 
produced by intense rf excitation pulses and stands ready 
to normally receive the signal. 
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II. 2 Measurement of Spin-Lattice Rela xation lime Tq 

1 28 

Measurement of T-^ "by the pulse technique s is des- 
cribed here briefly. 

SO 

Bloch's equation predicts that the growth of magne- 
tization along the Z-axis would be given by (see eq.(l.20) 

M z (t) = M q + [M z ( 0)-M q ] exp (-t/T-^ (2.2) 

where M Q is the equilibrium magnetization. From this equa- 
tion it is seen that a plot of log [M - I'i (t)] versus t 
gives a straight line with a negative slope equal to 1/T^. 
Therefore, experimentally, T, can be measured by first 
disturbing the system from equilibrium and then monitoring 
the growth of the Z-component of the magnetization as a 
function of time. [ M Q -M g ( t ) ] can be measured by applying 
a proper combination of % and tc/ 2 pulses. Two different 
pulse sequences were used, for convenience, in the present 
work to measure short and long T^. The procedure is 
described below: 

( i ) Measuremen t of Short Ti ( T-^ < 0.1 sec, as for B^-Ar, 

H 2 ~W 2 and E 2 -K" 2 0 mixtures) 

In this case a x-t-Tc/2-t ' -% pulse sequence was used. 

The first m pulse tipped the magnetization M to the 

z 

negative Z direction and its growth from -M 0 to +M q along 


the Z axis was monitored by the -jc / 2— t ’ — pulse sequence 
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applied at a time t after the first % purse. t* was kept 

constant G* o .5 msec) throughout the experiment so that the 

amplitude of the echo is proportional to the value of M 

preceding the %/2 pulse whereas t was varied continuously. 

M for t > 10 T, was taken to be the value of M . The 
z 1 o 

tc— t — tt/2— t 1 -it pulse sequence was repeated after every time 
interval T . which was taken to be about 10T, (T 10T-, ) 

JL 0 P b TQ P b A. 

to ensure that the spins relax and attain equilibrium in 

between the successive repetitions and after every repetion, 

i.e. for a new value of t, M z (t) starts growing from -M q . 

The time t between the first . % pulse and %/2 pulse was 

measured with a Beckmann Model 7370 counter. In order to 

improve the signal to noise ratio a Boxcar integrator 

Model ho. CW1 of Princeton Applied Research was used to 

sample the detected output from the receiver. The Boxcar 

was triggered by the n/2 pulse and the sampling gate width 

of the Boxcar was adjusted to sample the echo height. To 

avoid any possible integration distortion effects of the 

Boxcar integrator, the time t between the first n pulse and 

%/2 pulse was swept very slowly and the integration time 

constant RC of the v Boxcar integrator was chosen to satisfy 

31 

the following relation between the repetition rate, T, 
of the pulse sequence and the sampling gate width, t , of 
the Boxcar 

I- RC « X T, 

°s 1 


(2.3) 



18 


where a is the inverse of the rate at which the separation 

t between the first % pulse and %/2 pulse was varied and 

T^ is the spin-lattice relaxation time of the sample used. 

The output of the Boxcar, which gives the magnitude of the 

echo and hence the recovery of M , was plotted on a Bausch 

and Bomb strip chart recorder Model V0M6. The M z values at 

different times were identified by using an event marker 

on the recorder. T, was obtained by plotting log[M -M ( t) ] 

1 u z 

versus t. 

Here it should be noted that if the interval t between 
the first % pulse and n/2 pulse is varied not continuously 
but in steps manually, the relation given by eq. (2.3) need 
not be satisfied. But in that case one has to wait for a 
long time before changing over to a new t till a large 
number of repetitions with the same interval between the 
first n and m/2 pulses enter the Boxcar and the Boxcar shows 
a steady output. 

( i i ) Measurement of long relaxation t im es Ti 
(T T > O.l sec as for CHlcTJ' 

If T^ is measured according to the scheme written in 

(i) above, one has to wait for a long time of the order of 

lOT-^ between any two successive tt— t— tc/ 2— t 1 —vc pulse sequences. 

Por long T-^ it is a time consuming process and, to save 

32 

experimentation time, the following technique was used to 
measure T^: 
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Here the pulse sequence [ tc/2 burst] -t-Tc/2-t ' -tc was 

used. The tc/2 pulses in the hurst brought the magnetization 

to 0 and the growth of M from 0 to +I-i o was monitored by 

the Tc/2-t'-Tc pulse sequence. The height of the echo following 

the ix pulse was measured and it gave the value of F/ ( t) just 

before the tc/2 pulse. To ensure that M really started 

z 

growing from zero a large number (^lOO) of tt/2 pulses (a tc/2 
burst) were applied in a time (duration of the burst) much 
smaller than or . t' was kept constant and t was 
varied from 0 to about 10T-, . The value of M (t) for t ^ 1011 

X Z X 

was taken to be the value of M . It should be noted that the 
spin system need not go to equilibrium before each pulse 
sequence. Since the 7x/2-burst serves to ensure that no Z 
component of magnetization remains , it can be applied imme- 
diately after the echo and thus the [tc/2 burst ) -tc/2 --tc sequence 
is repeated after a time interval much smaller than 10 T^. 

The measurement of time t and recording of the signal 
was done exactly in the same way as described in (i) for the 
measurement of short The plot of log [M Q -M (t)] versus t 

is again a straight line with a negative slope equal to 1/T^. 

1 1 . 3 T he Spin-Ec h o Apparatus 
( i) General remark s 

The spin-echo spectrometer built by Raj an et al."^’^ 

QO 

based on the standard techniques as described by Hahn and 
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Clark was used in the present work. It lias an operating 
frequency of 30 MEz and makes use of phase sensitive detec- 
tion. The timing circuit of the spectrometer which usually 
gave n-ft/2~n pulse sequence was* however, modified to give 
the pulse sequence [n/2 burst] -% /2~% also. Thus, the 
spectrometer could be used to measure short as well as long 
relaxation times. Circuit details are given in Appendix. 

A block diagram of the spectrometer is shown in Pig. 2. 1'he 
major items are described below briefly. 

( i i ) T iming Circuit and the pulse sequ e ncer 

( a) 7i— ix / 2— it Pul se - sequence 

This pulse sequence was generated by using a combina- 
tion of the following three units: (i) a waveform generator 
lek 162 (shown as unit 1 in the Pig. 3), (ii) a pulse 
generator Tek 163 (unit 2), and (iii) another waveform 
generator Tek 162 (called as slow-sawtooth generator and 
shown as unit 3 in the Pig. 3) which is similar to the wave- 
form generator stated above in item (i) but modified to give 
a slow running sawtooth. A block diagram of the pulse 
programmer is shown in Pig. 3. 

The waveform generator (unit 1) operated in recurrent 
•mode produced (i) a triggering pulse, at a changeable repeti- 
tion rate, which was used to produce the first % pulse and 
(ii) a sawtooth. The sawtooth started running down at the 
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same time when the pulse began and had a width equal to the 
interval between the triggering pulses, i'his sawtooth was 
used to trigger the pulse generator (unit 2) which* within 
the duration of the sawtooth* gave a pulse of adjustable 
width and height at a time depending upon its triggering- 
level which can be adjusted manually as well as by the help 
of an external sawtooth. The slow-running sawtooth from 
the slow-sawtooth generator (unit 3) was used to serve this 
purpose. By using the slow-sawtooth it was possible to 
change the triggering level and hence ‘die time between the 
first Tt pulse and the pulse from the pnlse generator (unit 2) 
continuously. A reset mechanism was incorporated in the 
slow-sawtooth generator to enable one to stop and start the 
rundown at any time instead of waiting for the rundown to 
be completed of its own. 

The pulse from the pulse generator (unit 2) was t ak en 
to a pulse sequencer (Fig. Al) to produce two pulses ifhose 
separation was controlled by the pulse width of the input 
pulse. The pulse from the waveform generator (unit 1) and 
the second pulse from the pulse sequencer were added 
together in a mixer circuit. These tiro pulses were, then, 
taken to the 180° channel of the pulse width generator whereas 
the first pulse from the pulse sequencer went to the 90° 
channel of the pulse width generator. 
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( "b) |~7 e/2 burstl-Tc/2 -jz Pulse sequence 

In this case a collection of Tektronix waveform 
generators Tek 162 and pulse generators Tek 161 and lek 163 
were used. A block diagram of the pulse programmer is shown 
in Pig. 4- Positive pulses from the waveform generator 
Tek 162 (unit 1 in Pig. 4) operated in recurrent mode 
trigger the pulse generator Tek 163 (unit 2) which gives 
pulses of adjustable height and width. These in turn 
trigger another waveform generator Tek 162 (unit 3) to give 
sawtooths of changeable duration that are fed to a pulse 
width generator (unit 4). Its delayed pulse output was 
used to gate another waveform generator (unit 5) which 
gave a burst of very short pulses during the gating period. 

The number of pulses in the burst can be adjusted by chang- 
ing the duration of the gate as well as the repetition rate 
of this unit. The pulse from the pulse generator (unit 2) 
was taken to pulse sequencer (Pig. Al) to produce two pulses 
whose separation was controlled by the pulse width of the 
input pulse. The burst of short pulses from the last waveform 
generator (unit 5) and the first pulse from the pulse sequencer 
were added together in a mixer circuit (Pig. 4) and then 
given to the 90° channel of the pulse width generator while 
the second pulse from the pulse sequencer was given to the 
180° channel. The time t between the ii;/2-burst and the n/2 
pulse was changed by changing the repetition rate of the 
unit 1. 
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Block diagram: pulse programmer for 
producing TT-n/?-n pulse sequence- 











r ig. 4 Block diagram: pulse programmer for producing 
[rr/2 burst] - tt/ 2 -rr pulse sequence. 













26 


The two 90° and 180° channels of the pulse width 
generator (Pig. A2) were identical and were designed using 
phantastron circuits. The range of the pulse width could 
be varied from 0-70 psecs continuously ( -which can be further 
increased by changing EC appropriately) . These pulses were 
amplified in a pulse amplifier (Pig. A3) and then sent to 
the transmitter to gate it. 

( iii) Transmitt er 

The transmitter was designed to deliver high pulsed 
power output at 30 MHz with short rf rise ana fall times. 
Class C operation was preferred in certain stages to meet 
the power requirements. A crystal controlled oscillator 
(Pig. A4) at 10 MHz, housed in a rf leak tight copper can, 
was used as a source generator. The 10 MHz source signal 
was gated by the pulses from the pulse amplifier using a 
triode 7077 and then amplified (Pig. A3). The amplified 
output was trippled (Pig. A6) and then further amplified by 
a push-pull stage. The output from the push-pull stage 
was fed to the final power stage at 30 MHz (Pig. A7) . The 
power amplifier produced reasonably rectangular rf pulses 
of 1200 volts peak-to-peak amplitude. 

To make use of the phase sensitive detection the 
following arrangement was done. The output from the 
master oscillator was passed through a phase-shifter Ad-Yu 
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type 559 hi which was a commercial delay line and produced 
a delay of 0.4 psecs. The delayed output was passed through 
a trippler (Pig. AO) and then coupled to the last stage of 
the high gain L.E.L. amplifier in the receiver part of the 
spectrometer for phase sensitive detection. 

( iv) Receiver 

The voltage level of the signal inauced in the sample 

coil was first raised by using a preamplifier (Pig. A9) 

3 3 

which was similar to the design by Clark and 'makes use of 
tube 7722. The voltage gain of the preamplifier was about 
10 and the output impedance about 80 ohms. The signal was 
further amplified by an RCA 3005 integrated rf amplifier 
having a power gain of 25 db and output and input impe nances 
of 50 ohm. The circuit of the integrated amplifier is 
shown in Pig. A10. The design was taken from the RCA appli- 
cation note ICAKf-5022. The amplified NMR signal was ? then* 
fed to the receiver which was a L.E.l. amplifier of the 
Varian Associates s model Ho. 1MN2-30-15-50, having a centre 
frequency of 50 MHz s 3 db’ bandwidth of 10 MHz and input 
impedance of 50 ohms. One of the two output terminals 
provided in the L.E.L. amplifier lies just after the 
detection and has a gain of 90 db maximum while the other 
one is after a video stage and has a gain of 115 db minimum 
at 0 volt bias. 
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The L.Ei.L. amplifier used a 1H295A diode detector. 
Ordinary diode detection is not linear over the entire region 
and the non-linearity becomes important while dealing with 
small signals. However, a diode detector can be made to 
operate in the linear region by using phase sensitive detec- 
tion. In phase sensitive detection the diode detects the sum 
of reference voltage that biases it and the signal. So 
using the phase sensitive detection the diode detector of 
the L.E.L. amplifier was made to operate in the linear region 
by choosing the reference voltage properly so that the operat- 
ing point fell in the linear region. To find the proper 
biasing (reference) voltage for the diode an experiment was 
performed, for a given setting of the spectrometer, the 
reference voltage was varied and the output signal level was 
recorded each time. The output signal strength was plotted 
as a function of reference voltage (Pig. 5) . Prom the Pig. 5 
it is seen that the reference voltage biasing the diode 
detector must be kept between -1.5 and -2.5 volts for 
linear detection. To avoide distortion the signal input to 
the diode detector was kept below 1/10 of the reference 
vol tage . 

( v) Sample coil, sample holder and the gas handling system 

The spectrometer makes use of a single coil which 
optimizes the transmitting and receiving modes. The schematic 
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of the spectrometer containing the sample coil is shown in 
Fig. 6. Depending upon the two different types of measure- 
ment in the present work, two separate sample coils were 

used; one for OH- Cl where the measurements were confined in 

o 

the low density region (0.2 < P <6.5 amagats) at room 
temperature and the other for T-^ measurements in the hydrogen- 
mixtures in which case the densities of the mixtures went 
as high as 70 amagats and the temperature varied from 300K 
to 600E. In order to improve the signal to noise ratio for 
weak signals in the low density region a sample coil of 
bigger dimensions was made. The sample holder carrying 
this bigger coil did not leave any space between the pole 
pieces of the magnet to incorporate a heater. So, for high 
temperature work (EL^ mixtures), a smaller sample holder 
was used. The details of the coils and sample holders are 
given below. 

(a) For measurements at low pressures ( CH-^ Cl as sample ) 

The coil used in this case was of length 3.5" and 
diameter 0.8" and consisted of 14 turns of 18 SWG copper 
wire. The enamel insulation was completely stripped off 
and the wire was thoroughly cleaned before use. The coil 
was housed in a sample holder (Fig. 7) made of a thick — 
walled Be-Cu vessel (1,5" o.d., 1" i.d.) with brass plug A 
sealed to the vessel with a teflon washer. A glass tube 
which could fit between the coil and the wall of the sample 
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Schematic : coupling of RF head. 
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holder was used to insulate the coil .from the sample holder. 
After fixing the return end of the sample coil to the plug A 
with a screw, the gLass tube was put in position such that 
the other end of the coil came out through the small central 
hole of the tube as shown in fig. 7. ihis end was taken out 
through a 1/8" bore in the plug B. Thin glass capilliary 
tubes were used to insulate the wire from the body of the 
plug. The plug B was then tightened to make pressure seal. 

The rf head was sealed against pressure using teflon washer. 

The lead of the coil coming out of the rf head was soldered 
to a BBC connector. 

(b) for measur em ents at high pressures ( EU mixtures as sample ) 

17 

A sample holder made by Rajan et al. was used in this 
work. The coil was of length 3.8 cm and 0.8 cm diameter and 
consisted of 15-16 turns of a cleaned 20 SWG copper wire. It 
was housed in a thick-walled Be-Cu vessel (o.d. 2.7 cm and 
i.d. 1 cm) with plugs A and B (fig. 8) sealed to the vessel 
by using the technique of AE cone fitting. first the return 
lead of the sample coil was screwed to 'the plug B. Then B 
was sealed in its position by tightening the nut E. The 
other end of the coil was taken out through a bore in A and 
A was then tightened to make a pressure seal. The rf head 
was sealed against pressure at room temperature using teflon 
washer. Glass tubes were used to insulate the coil from the 
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Fig . 7 


Sample holder used for CH 3 CI . 
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sample holder exactly in the same manner as described in 
the preceding paragraph in (a). 

Pressures above 4 PSIG were measured with Bourdon 
pressure gauges which were calibrated using a dead weight 
gauge and below 4 PSIG were measured with a mercury manometer. 
In the experiments with CH^Cl the use of rubber was totally 
avoided since CII^Cl reacts with rubber. 

Ihe schematic of the gas handling system is shown in 

Pig. 9. 

( vi i ) Temper at ur e_ measurements 
( a) Above ro om temperature 

High temperature work was done only with the smaller 
sample holder shown in Pig. 8. A commercial 1000 watts non- 
magnetic spiral heater element was wound on its outside 
non-induct ivel y . The heater was insulated from the body of 
the sample holder by a thin mica sheet. The windings of the 
heater were held in place by using a liquid porcelain 
called ' sauereisen 1 . D.C. was used to heat the sample holder 
since small 50Hz noise was observed when A. C. was used. The 
high-pressure vessel along with the heater was kept inside 
a vacuum jacket to keep the convection and conduction current 
losses to a minimum. In order to maintain the sample holder 
at 600K about 120 watts of electrical power had to be 
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dissipated. The leads from the heater (and the thermocouple 
also) were taken outside the vacuum jacket through Kovar 
seals . 

A Chr o me 1 -Alume 1 thermocouple was fixed 
away from the heater on the top plug of the sample holder to 
measure the temperature. The temperature measurement was done 
with a Honeywell potentiometer model 2705B which gave an 
accuracy of +1°, Ho special effort was made to control the 
temperature automatically because once the current was 
adjusted the temperature was steady within a degree after 
reaching equilibrium. However, at least 8~10 hours were 
allowed before starting the T^ measurements to ensure that 
the sample was at equilibrium with the sample holder. 

( b) Below room temperature 

The temperatures in the range 300-20 OK were produced 
by using a brine, spirit+liquid nitrogen, the temperature of 
which depends upon its composition. The more is the liquid 
nitrogen, the lower is the temperature. Hence by varying 
the composition a. range of temperature from room temperature 
to the freezing point of spirit can be achieved. The sample 
holder was put in a Dewar flask filled with the brine and 
the flask was covered with glass-wool. ¥hen the system had 
attained equilibrium a temperature variation of only 0.5 
degrees per hour was observed. The temperature was measured 
by the Chromel-Alumel thermocouple as described in the 
preceding paragraph. 
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( vi ii ) Preparation of sample 
( a) Purity of g ases 

H 2 , N 2 and used in the present work were of research 

grade and were obtained from Matheson G-as Co. Argon was 
obtained from Indian Oxygen Ltd. and was of 99. 99 percent 
purity. Methyl chloride gas was supplied by Mettur Chemicals, 
India and was of 99-95 percent purity. However, it was 
suspected from the measurements that CH^Cl contained some 
impurities. So it was further purified in the following way; 

CH^Cl gas cylinder was connected to a 20 ft. long y n 

51 

i.d. copper tubing packed with Zeolite molecular sieve 

( 3A, appropriate for absorbing oxygen and water vapour) and 

which was made into a coil of 1 ft. diameter. The other end 

of the coil was connected to a stainless steel container 

through pressure valves. The copper tubing and the container 

were thoroughly evacuated. CH^Cl gas was, then, passed 

through the copper tubing very slowly and stored in the 
oo 

continer. The molecular sieves were Regenerated in the 
beginning as well as intermittantly after some time intervals 
much shorter than the saturation time of the sieves by 
heating the coil to 150°C and then cooling it. During 
the process of regeneration, gas flow was stopped and the 
coil was continuously evacuated. The oxygen free gas, 
stored in the container, was further purified by freezing it. 
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To do this the stainless steel container was put in a dewar 
containing ethyl alcohol + liquid nitrogen mixture. The 
composition of the mixture was adjusted to maintain its 
temperature around -40°C, a temperature he low the boiling 
point of CH^Cl. The residual gases in the container were 
pumped out. The container was further cooled to -132 C, 
a temperature below the melting point of CH^Cl by adding more 
liquid nitrogen to the cooling mixture simultaneously pumping 
out the residual gases from the container. The container 
having a high purity solid CH^Cl was then brought to room 
temperature . 

( b) Mixing of gases and composition d et exam ination 

A mixture Hp-X, where X is Ar, Np or NpO, was prepared 
in a separate container. First the container was filled with 
the low concentration component to a pressure which would 
give approximately the required composition at the total 
pressure of the container. Then the high concentration 
component was introduced to bring the total pressure to the 
desired value. The composition of the mixture was estimated 
from the partial pressures and was checked by comparing the 
signal strength (echo-height) of the mixture with that of 
pure Hp under identical conditions of pressure, temperature 
and spectrometer adjustments. 

The. signal strengths for pure Hp and Hp-X mixture were 
corrected for Tp effects. The height of the spin-echo A (2t*)» 
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where t' is the time between the n/2 and. % pulses, is given 

A(2t') = A(0) exp [- “"] (2.4) 

2 

where diffusion effects have been neglected. In terms of 
the amagat density p where an amagat is 2.63 x 10^ molecules/ 
cc, A(0) can be bitten as, 

A(0) = CP • (2.5) 

with C as a constant depending upon the number of spins per 
molecule and temperature. At high densities it can be 
shown that 

3-l/P = ? 2 /p (2.6) 

and eq.. (2.4) re duces to 

A(2t») = C P exp (-2t'/T 1 ) (2.7) 

or CP = A( 2t 1 ) exp ( 2t ' /\) (2.8) 

( 0 3 ) I T . 

Thus the quantity — 7 gives the fractional concen- 

1 0 p; Pure H 2 

tration of in the mixture. 

It should be noted, as is clear from the eq. (2.7), that 
for the samples where >> • 2t 1 , A( 2t ' ) % C P but for the 
cases such as H 2 where a, 2f, the exponential factor 
called the ' I 2 correction' is important. 



CHAPTER III 


PROTOH SPOT-LATTICE RELAXATION IE 
H 2 -Ar, H 2 -N 2 AND H 2 -IT 2 0 KEXOTRES 

III.l General Remarks 

The interhal degrees of freedom of a H 2 molecule corres- 
pond to (i) the electronic (orbital and spin), (if) vibra- 
tional, (iii) rotational and (iv) the nuclear spin states. 
Since only the ground electronic and vibrational states are 
occupied up to a temperature 1000K, the molecule can be 
considered, below this temperature, as a rigid rotor having 
an additional degree of freedom corresponding to the nuclear 
spin states. Thus, the state of H 2 molecule can be described 
by the rotational quantum number J, nij and nuclear spin 
quantum numbers I, m^.. 

Further, from the requirement that the total wave- 
function of the hydrogen molecule be antisymmetric with 
respect to the interchange of the two protons, it is 
concluded that, in general, hydrogen gas consists of txro 
types of molecules, one having the spins of the two nuclei 
parallel (total spin 1=1) and occupying the odd rotational 
states (J = 1, 3, 5 etc.) and the other having the nuclear 
spins antiparallel (total spin 1=0) and occupying the even 
rotational states (J = 0, 2, 4 etc.). These two are called 
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ortho and. para modifications respectively. The number of 
molecules occupying a J state is given by the Boltzmann 
distribution 

Pj = gj exp ( -Ej/kf) (3.1) 

where Pj is the fraction of the total number of the molecules 
in the gas occupying the J state, k is the Boltzmann constant 
and I is the absolute temperature, gj is the degeneracy of 
the states 

gj = ( 21+1) ( 2 J-i-l) = 3(2J+l) for 1=1, i.e. J odd (ortho-H^) 

= (2J11) for 1=0, i.e. J even (para-TL,) 

(3.2) 

Ej = = J(J+l) k0 R (3.3) 

and is the energy of the J state. Here I = moment of inertia 
of the Hg molecule, -fx = Plank's constant. 

e R = 2 = 85.3°K for H 2 (3.4) 

The values of Pj for different temperatures are given in 
table 3.1. 

The ratio of ortho- to para-molecules in a sample of 
hydrogen at thermal equilibrium is given by 

I 3(2J+1) exp [-J( J+l) 0 p/T] 

J=odd 

'l (2J+1) exp [-J(J+1) 0 p/T] 

J=even 


r 


(3.5) 



Table 3.1. Fractional population of the rotational states for 
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From the above expression it is seen that at temperatures 
l 1 >> e R , r_*3, i.e. about 75^ of the total number of molecules 
in the sample are ortho-molecules and the rest are pai’a- 
molecules. As T + 0, r +0, i.e. as the temperature goes down 
the ortho molecules get converted to para form. But this 
conversion does not occur so quickly because the probability 
of transition of ortho-E^ to para-E^ which involves change 
of total nuclear spin is very small and it takes a time of 
the order of years for the equilibrium to be established 
between the ortho- and para-forms. 'therefore, normal H 2 
gas is considered as a mixture of two separate gases, nameljy 
ortho-H 2 and para-B^. A mixture H 2 ~X, thus, should be 
treated as a mixture of three species, viz. ortho-H 2 , para- 
H 2 and X. 

If Ptt and p v is the number of and X molecules per 

I±2 4 . 

unit volume respectively and and represent the number 
densities of ortho- and para-H 2 then the density, p^,, of the 
mixture is given by 


P T “ P H. 


+ p 


X 


where 


H, 


p 0 + p P 


(3.6) 

(3.7) 


the fractional densities of the components in the mixture 
can be defined as 


p x 


X- 


- i brr 
P T H 2 


PH. 

4 

P/p 


PO 


p 0 

Pj 


(3.8) 



( 3 . 9 ) 


and for normal Hg 

p 0 = (3/4) P H 

Only tlie ortho-Hg molecules give NHR signal since the 
para-Hg does not have a nuclear spin moment. A sample 
containing N number of molecules gives a signal having the 
strength proportional to ^NI(I-KL) = . 

I I I . 2 Review of the Theory 

The theory of nuclear spin relaxation in and its 
mixtures is- given in detail elsewhere'*’ 5 '^ and references 
cited therein. Bloom and Oppenheim have obtained expressions 
for the relaxation rates explicitly in terms of the inter- 
molecular interactions by expressing the correlation func- 
tion of the intramolecular interactions in terms of the 
correlation functions of the intermolecular interactions. 

Here, a brief review of Bloom-Oppenheim theory"^ for I s 
presented and expressions for proton relaxation rates in 
mixtures useful in the analysis of the present data are 
given . 

C orrelation functions of the intramolecular interactions 

In a sample containing N molecules of hydrogen gas kept 
in a static magnetic field H Q along Z direction,, the fluc- 
tuating spin-rotational and dipolar intramolecular interactions 
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produce a small time dependent perturbing term sfC in the 
total Hamiltonian which is given by"*"’"*"^ 



N 



Yffl 1 

2 




( 6 ^)* 


yhH" 





( 3 . 10 ) 


where the index i denotes the ith molecule and 

y = 2.69 x 10^ sec - ^ gauss'”"*" is the proton gyromagnetic 
ratio 

H’ = 26.75 gauss is the coupling constant for intramole- 
cular spin-rotational interaction 
H" = 53.86 gauss is the coupling constant for the intra- 
molecular dipolar interaction 

and 1„ and S«,™ (& =1,2) are the functions of molecular 
&m 

rotational operators and nuclear spin operators respectively. 


L 10 

= V"2 j 

z 

S 10 

= V2 I z 


1+1 

II 

II 

J + i J 

3C y 

s i+i 

= *+ = *x ± 

i3 y 

0 

CM 

O 

CM 

il 

(n.) 

S 20 

= (2/3) ^[31 

2 - 1 ( 1 + 1 )] 

J 2+l 

— Y 

2+1 

(») 

S 2+l 

= II, + I, 

Z + + 

— “ T"" 

J z 

'2+2 

= Y 

2+2 

(n) 

S 2+2 

= I 2 



where Y 2 m (ft), m = 0, +1, +2 are the spherical harmonics of 
order 2 of the orientation q of the molecular axis. 
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By using the general theory'*' the proton spin -lattice 
relaxation time can be written as 1 ’ 10 ’^ 

= I - 2 |- 2 J 11 (a) I ) + (i|S) y 2 H" 2 + 4J 22 (2 Wi )] 

(3.11) 

where 

J i.m = f e ' iUt (3 ’ 12) 

CO 

is the Fourier transform of G-^t) and is known as the spec- 
tral density, and 

4 m (7 + 4m ( °b (3.13) 

is the correlation function of The bar represents an 

ensemble average. It was well established by Schwinger , 

4 1 

Needier and Opechowsl-ci , Abragam and was also shown experi- 

37 

mentally by Hardy that the intramolecular correlation func- 
tions can be expressed as exponential functions of time which 
decay with a time constant t known as correlation time of 
the intramolecular interaction. Johnson and Waugh pointed 
out that out of different possible correlation functions only 
two, G 1q and G 2q , are important. In the short correlation 

time limit since the values of the spectral densities of G ^ 

2,0 

evaluated at frequencies or are the same as those 
evaluated at zero frequency the eq. (3.11) reduces to 

?jr- = -b 2 H' 2 J lo (0) + 12 tc y 2 H" 2 J 2 o (0) 


(3.14) 
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The correlation functions G 0 can he written in terms 

of the matrix elements of 1. between the rotational states , 

20 

in the folloxriLng way 




= ji I vCJH,t/J o « o ,0)(JM|I )lo | JM) 

^(3 x o 0 d 


= i._ P .T M [ f n (Jn)] 2 C(J 0 £J 0 ’ M o 0 ) I w( JM| t / J n M n j 0 ) X 


J 0 M 0 ‘W 0 


ll 


O 0 * 


[f £ (j)] 2 C(J£J 5 M0) (3.15) 


where 

( JH[L Jl0 j JM) = [f £ ( J# C ( J£ J 5 MO ) 
is the matrix element of l £o? C( J£J; MO) is a Clebsch-Gordon 

'TjO 

coefficient (Rose^ ) 


%( J) 


f 2 ( J) 


2J( J+l) 


j(j+d 


4n ( 2 J-l) ( 2 J+3) 


3114 Pj o H o ' 2 J o«- Pj o ~ 


exp (-Ej /kT) 
YT2J+1) exp T-B jTkTT 


(3-16) 


(3.17) 


J - 


is the probability that the molecule is in the state J qS M q 

with energy Ej while Pj is the probability that the molecule 
d o °o 

has the rotational energy Ej • w( JM #t/J Q M 0 jO) is conditional 

d o 

probability that the molecule is in the rotational state J,M 
at a time t given that it is in the state J 0 sM 0 at t=0. 

^16 3X1 ^ ^2o are ca l cu -l a l e( ^ 6y assuming that the condi- 
tional probability w(JM$t/J 0 M jO) satisfies a temporally 
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homogeneous master equation 

If = w(JMjt/J o H o ,0)= ^ t A(JI4} J'M') w( J'M' ;t/J o M o jO) (9.18) 

where A( Jl; J'M' ) , J, M ^ is the transition rate for the 

transition J* ,M' J,M that is caused because of the inter- 
molecular anisotropic forces acting during collisions. 

Por J=J ! , M=M' 

A(J'M' 5 J'M' ) =- l l ) - l u J’ K, J'M’) (3.19) 

J^J' M M^M* 

Prom eq. (3.18) it is seen that a tensor polarization of order 
l in the state J for molecules initially in the state JMi 

O o 

and defined as 

m £ TTM = l C( JS,J }M0) w ( JM ; t / J M j 0 ) (3.20) 

^0 0 M 0 

satisfies the following relation 

m £ TT m (t) = I w( J' M' s t / J M ;0) l C(JAJjMO) A(JM,-J'M') 

o n o J'M' ° 0 M 

(3.21) 

The number of coupled equations* eq. (3.18)* is equal to the 
number of rotational states appreciably populated at the 
temperature of interest. 

The transition rates A(JM$J'M’) are calculated using 
'weak collision approximation' which means that the anisotropic 
in term'd ecular potentials are assumed to be weak in the sense 
that first order perturbation theory may be used to calculate 
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the transition rates A(JMjJ , M l ). If each term of the inter- 
molecular anisotropic interaction potential be expressed as 
a product of a lattice operator involving the intermolecular 
separation and a rotational operator involving the orientations;, 
of the molecular symmetry axes, which can be represented by 
spherical harmonics Y ( 0) of order X , then the transition 
probability A x (JMjJ’M' ) which contains the squares of the 
matrix elements of the operators between the states J' ,M' 
and J, M can be written as 

A x (JM,-J’M>) = Q x [C( J'X (3.22) 

where Q x (J, J') depends upon the radial dependence of the 
anisotropic intermolecular potential and involves the corre- 
lation function of the intermolecular interactions. Using 

38 

eq. (3.22) and the eqs. (6.23a) and (6.24) of Rose ", it 
can be shown that 

^ C(J£JjM0)A x (JMjJ'M*) = C(J , £J’ |M , 0)B £X (J, J’ ) (3.23) 

where 

(|J1+1)B^ ( J,J') = ^(J,J 1 )(“1) X+; '" J “ J ’[(2J+1)(2J , +1)]^ 

x W(J, J’ , J, J' j\£) -6 jj, l (|~~)Q X ( J M ,J’) 

(3.24) 

7n 

where W( J, J' , J, J'jXA) is a Racah coefficient (Rose^ ) . From 
these expressions and eq. (3.21), it is seen that the tensor 
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polarization nf TT satisfies the following relation in tlie 
^ J o <d 

weak collision approximation 


m JJ M ~ I 2 0 (J, J' ) 
oo J 


Jt, 


• ' mj,J o M o 


(t) 


(3.25) 


where 


B, 


= l E (J,J*) 


( 3 . 26 ) 


Solving the eq. (3.25) with the initial conditions 


m JJ M 
o o 


C(J o £<V M o 0) 6 JJ 


and w( M >0) 


6 jj 0 6 mm 0 


and using eq. (3.15) the correlation functions G (t) are 

X* ^ 

obtained in terms of Q^(J, J')» i.e. in terms of the correla- 


tion functions of the intermole cular interaction, 


Correlation functions of the intermoleculax° interactions 


i'he inter mole cular potential between two molecules 
consists of isotropic and anisotropic parts. The important 
terms of the anisotropic interaction potential, v^, between 
two molecules labelled 1 and 2 can be written as a sum of 
two terms one depending upon the orientation of the first 
molecule and the other depending upon the orientation and the 
relative position of the second molecule also. If r be the 
vector distance of the centre of mass of 2 from that of the 

molecule 1 and ej the angle which the synfn$ t£y jff 

Ci'P' JN'vl. 


n> 


7 62235 
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(3.27) 

(3.28) 

(3.29) 


in which. P^cos 6'^) is the Legendre polynomial of order 2 of 

the angle 9^, "^q^l^ and ^2q^P are spherical harmonics of 
? * 

order 2 and ^ and 9 .^, are , respectively;, the orientations 
of the symmetry ares of molecules 1 and 2 with respect to r. 
b^"0(r) fc(2)^ r ) are function of r and a^ are constants. 

For quadrupole -quadrupole interaction a 0 =6, a 1 =a_2 = "4? 
a 2 =a_ 2 =l and 



4n:Q 1 Q2 



(3.30) 


where and are the quadrupole moments of the molecules 
1 and 2 respectively. The potential written in eq. (3.27) 
is the most general one if both the molecules are ortho-Hg 
molecules in their lowest rotational state. If the molecule 2 
is a para-H 2 or inert gas molecule, v^ 1 ^ represents the 
interaction potential. If the molecule 1 is an ortho-E^ 
molecule in any J state and the molecule 2 is another molecule 
having permanent quadrupole moment, such as v , as written 

in eq. (3.27), gives the leading terms of the anisotropic 
intermolecular potential. 
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In the laboratory reference frame, in which tlie Z axis 
is along H Q , v^'^ are written as" 1 ' 2 ’ 


v A (1) =|S b^p) I _ Y^fap T^fe) 


P=-2 


2p3“-l' 2p, 


(3.31) 


v 


A 


(2) = (||2)i b (2 ) (r) . j ( -1) ^C( 224*,m, m 1 , p) Y pm (P ) Y ?rn , >>) • 


m m' 
P 


"W 0 * 


2m 1' 2m* c 

(3.32) 


where P-^, ® are the orientations of the molecules 

1, 2 and the vector r , respectively, in the laboratory frame. 

Since v^"^ involves rotational operators of the molecule 1 
( 2 ) 

only and v^ v ' involves rotational operators for both the 
molecules, it is seen that there is no interference between 
the contributions, from and v , to Q x ( J, J 1 ) , so that 

one can write 


<^(J,J') = Q^^U, J') + q/ 2) (J, J'),X —2 for H £ (3.33) 

where J' ) and ^ 2 ^ ( J, J' ) arise from v^ 1 ^ and v^ 2 ^ 

respectively. The effect of v^^ is to produce transitions 

from the state J* ,M' to J,M in the molecule 1 and that of 

( 2) 

v^ v is to produce transitions in the molecule 1 from J' ,M ! 
to J,M while the molecule 2 simultaneously undergoes transi- 
tion from the state J ,H ,M' n to J", M" . The transition rates 
A^^CJM, J'M') and A^ 2 ^ ( JM, J* M' ) , for the molecule 1, produced 
by v^ 1 ) and v^ 2 ^ respectively can be written using the first 
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order perturbation theory as 

» ± w t T , t 


( JM§ J' M' ) = ~ \e JJ< ( J' M' | vT^ ( o ) | JM) ( JM | 2 ^ (t 7 jj*ip) dt 

xl «— m *”* 


A 


u'C 


(5.54) 


= 5” ( l^ ) t C(J ' 2J » OO )] 2 C c ( J,2J ^ M ^)] 2 0 (l) ( a, jj») 


where j (l) ( UjJ1 ) = f k (l) (t) dt 


(5.55) 


fan 


JJ' 


] J ~ E J< 


in which 


.( 1 ) 


k VJ - ; (t) = (H-l) b^(r(0)) Y (a (0)) b U; (r(t))T (ft(t)) 


TIT 


•2-n' 


(5.56) 


is the correlation function of b^(r) Y« (n ) 

2-p 


A^( JM ; J'M') = I I A[( JMjJ'M' )( J" r M ,H )l (5.57) 

jMjm M »*M •* * 


where 


A[ ( JMj J’ M' ) ( J"M n 5 J ,u M ,n ) ] = l (||2L) j (2) (w 

mm * ^ dd ;0 d 

p 

X [c'(224sm,m')] 2 ^rhr|(JM|Y 2 m (Si)| J ' M >)| 2 

x- | (J" M"|7 2m ,(n 2 ) 2 (3.38) 

and represents the rate of simultaneous transitions produced, 
in the molecules, by . 
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Here dloj 


JJ' f J"J 


n,T 1 1 1 “ (Bj— Ej,) (®.tu “ ■‘At*" ) 


J' 


(3.39) 


( 2 ) 

and y 1 ( ^ j„j„i ) is the Fourier transform of the correla- 

tion function k^ 2 ^(t) of b^ 2 ^ ( r) (ft) , where 


k (2) (t) =(U-1) b^"(r(0))T 4 (0 (0)) b (2 ^(r(t)) ? 4 (3(t)) 


(3.40) 


(a ) 

values of j Viy («) 5 q = 1,2, have been calculated by Bloom 
et al.^^ using Constant Acceleration Approximation ( CAA . 
Summations in eq. (3.37) are performed by the method given 
by Rose. ^ 


Using these expressions for A^^ ( JTi$ J’ M l ) and A^ 2 ^ ( JM$ 

J'M* ) and eq. (3.22) where X = 2 for E , ( JJ 1 ) and 

( 2 ) 

Q 2 M JJ') are obtained, in terms of the spectral densities 
of the correlation functions k^^(t) and k^ 2 ^(t) of the 
intermole cular interaction, as 


Q 2 (1) (J,J’) = (|^r ; )[C(J , 2J|00)j 2 o' (1) (“ JJI ) (3.41) 

and Q 2 (2) (J,J») = 5^CC(J2J‘ 5 00)] 2 ^ JIt( r j" ! 

x [ C( J ,M 2 J” j 00) ] 2 (3.42) 


Using these expressions and the relations 'that preceded G- 

Xf. O 

are obtained which with the help of eq. (3.14) yield an 
expression for 1-^ in terms of the spectral densities j^^(u>.jj 


) 
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( P) 

and 0 '(w JJt 
tions . 


jitjm) of the intermoleculau correlation func- 


In the- present work, mixtures of H 2 with Ar, 11 2 and h T 2 0 
have "been studied. H 2 '--Ar collisions cannot cause any transi- 
tions in the J states of ortho -Ii 2 "because Ar is a spherically 
symmetric atom. H 2 ~E" 2 and H 2 ~N 2 0 collisions also cannot 
produce any inter J transitions in the rotational states of 
ortho-H 2 "because of the energy mismatch between the J levels 
of H 2 and B" 2 (or B^O) for allowed transitions. Hence the 
case of infrequent transitions between J states of H 2 will be 
considered in detail and expression for will be given. 


Case of infrequen t transitions between states of different J 

For the case where the transitions between different 
rotational states occur much less frequently than those between 
the states with different M within the same J manifold, it 
can be shown that 



Q 2 (J, J' )<< Q 2 (J, J) j 5^ ■ J ! 

(3.43) 

and 

j (1) <“jjO = j (1) <°) Sjj. 

(3.44) 


J"0 = 3 (O (0) 6jj.Sj.ijim 

(3.45) 


£>qs, (3.41) and (3.42), then reduce to 


Q ( i ) f T T") 4'"^ . J( 1 ) _; ( l) ( f\\ 

Q 2 (t],d) _ 5 X2J^T)(2J+3) d ^ 


(3.46) 
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and 


q 2 (2) ( j, j) 


1 __ J(J+1) 

“ 25 % "(" 2 J-iT( 2 


C j (2) (0) 


(5.47) 


where 



J» (J"+l) 
T"2J ,, -]I)( 2 J h + 5j 


_ 5 / J( J+l) 

“ 2 \I 2 J— 1) ( 2 J+ 


/ 

3j)(5.48) 


and the prime indicates that the average is taken over the 
rotational states of the molecules 2 which collide with the 
ortho-H^ molecule (molecule l) . 

q io 

Bloom and Oppenheinr' ’ have evaluated expressions for 
0^(0) using constant Acceleration Approximation 1 ^ as 


j( < i)(o) = paliisaali i ( i) ( - q) 


(3.49) 


where p' is the number density of the gas, p, is the reduced 
mass of the colliding pair and a is a length parameter which 
is a measure of the range of the isotropic intermole cular 
potential. The index p is determined by the angular variation 
of the anisotropic interaction. When q=l, p=2 and when q=2, 
p=4 for H 2 . 


If the radial dependence of the anisotropic inter- 
molecular potentials v^^ (see eqs. ( 3. 27) -( 3. 30) ) could be 
written in the form 


,<n (ax) 


x 11 ' 


x' 


.n 


(3.50) 


where r = ax and x is a dimensionless variable, then one can 

Q 

write^ 
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Substituting eq. (3.49) into eqs. (3.46) and. (3.47), one 
obtains 

Q 2 ( 1 ) (J,J) = (^~ ) pa 4 ( 2 1 T( 3 | J )n (l) ( 2 ) [; 2 ^i 5 x' 2 T + Tr ] (3 ' 57 > 

and Q 2 ( 2 ) (J,J) =(-~~) p a 4 ( 2n[3u) *1 ( 2) ( 4) [fgjMj B J+5) ] C (3>58) 

Using eqs. (3-57), (3.58), (3.33) and (3.24), B £ (JJ) are 
obtained as 

t l(1)(2) + “r ® (2) (4)] 

(5.59) 

and B 2 (J,J) = ^ 2 ~ j ' 4i fr 2J+3) 7 ^ % ( J > J) ( 5.6° ) 


which, in turn, can be expressed in terms of quantities k Q and 
k^ defined as 


k-, 


6 tc 

5fi 2 

21 

50 nh‘ 


I (l) (2) 


( 2 ) 


(4) 


(3.61) 

(3.62) 


in the following manner 


and 


Bi( J, J) 
b 2 ( J, J) 


2 p a^( 2 ti Bu) 2 ( y. , an \ 
rrmuMT (k o ck i^ 


is given by eq. (3.60). 


(3.63) 


Using these expressions for B £ (J, J) and eqs. (3.15), (3.20), 
(3.21), ( 3. 22) -( 3. 26) and (3.43), and are obtained 
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whose Fourier transforms are given by 


10 


£o 
and. T, 


" a+l J p j f A ( J) l^T j7j}[ i » £ " 1 » 2 

is, then, written with the help of eq. (3.14) as 


(3.64) 


Pa^( 2 ( k : +Ck, ) -] 

o._ l_ = l y 2 H , 2 {j( J+l)(2J-l)(2J+3)} 

1 J 


(1 + 


1 


3H" 2 

H* 2 4 J 2 +4 J-7 


>] 


or 



1.02X10 13 <b( J) > 
pa^(27ipp) T ( k^+Ck^) 


(3.65) 

( 3 . 66 ) 


whe re 


h( J) 


J( J+l) ( 2 J— 1 ) ( 2 J+3) 
57.4 


(1 + 


4.74 


4 J 2 + 4 J-7 


(3.67) 


* 

and < > Represents an average over all the rotational states of 
the ortho-IL, molecule (molecule 1) . The numerical values of 
y, H' , and H" for have been used to find out the numerical 
factors in eqs. (3.66) and (3.67). 


E^-X mixtures (X = Ar, and N^O) 

The total interaction potential between ortho-E^ and X 
molecules consists of an isotropic part v Q (r) and an aniso- 
tropic part v^(r,e v )» where r is the in ter molecular separation 
and 0' represents the orientations of the molecular symmetry 
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axes with, respect to r (i.e. in body fixed frame). It is 
assumed that v (r) is given by Lennard -Jones potential 

v 0 (r) = 4e [(f) 12 - (f) 6 ] { 3.68) 


where s and a are the Lennar d Jones energy and distance 
parameters appropriate for B^-X pair. 


The anisotropic part of the interaction potential 
between and Ar molecules depends upon the orientation of 
molecule only , since Ar is a spherically symmetric atom, 
and hence can be written, in the form of v^ 1 ) as given in 
eq. (3.27). Bor E^-N^ and H^-N^O pairs the interaction 
potential v^ contains both v^"^ and v^^ because and 
K^O have permanent quadrupole moments. 


The long range attractive parts, v ^.( r,8 ' ) , of the 
anisotropic intermole cular potentials between H 2 and X can 
be written in terms of their polarizabilities, following 


Buckingham's results 


10 


as 


3U, U 9 a, 1 ( a -a ) , 

t v att^ r,e ' ^X=Ar = ~ ~2(Uj+U 2 ) a 2 ^ + 3 P 2^ C0S e i) 


^ v atfc^ r,e '^x=N, 


+ ] (3.69) 

T 


or NpO 


3U 1 U 2 r a i a 2 1 “^1 -r, , _ . V 

STU-j+Ug) ^ ^6” + 3 2 r 2 ( S 9 1^ 


! («„ -a x ) 

+ 3 a l E ~ P 2 ( cosQ ?) + ] 


I 


(3.70) 
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where only the leading terms have been retained and 8^ and 
are respectively the orientation of the molecular symmetry 
axes of molecules 1 and 2 with respect to 3*. and a ^ are 

the total polarizabilities of the molecules 1 and 2 while 
a ( j and are the polarizabilities parallel and perpendicular 
to the molecular axes . an ^ ^ ^ ave ^ lie un i"* :s °f energy 

and it was suggested by Buckingham^ that they might be taken 
as the first ionization potentials of the molecules. The 
first terms in eqs . (3.69) and (3.70) are the isotropic parts 
of the long range interaction and hence must be identical 
with the attractive part in eq. (3.68). Therefore it was 
chosen to fix 


5U 1 U 2 


a l a 2 


= 4e a 


3.71) 


It should be noted here that the values of U-^U^/^IJ^+TJ^) 
obtained for Hp-X pairs using the above equation are approxi- 
mately twice of those obtained following Buckingham's sugges- 
tion. This difference is due to the fact that the r - ^ term 


is only the leading term of the potential derived by Buckingham 
while in 12-6 potential it is the only term. 


It is assumed that the anisotropic interaction potentials 
. -12 

contain repulsive terms also which vary as r and have the 
same angle dependence as that of the attractive parts. Thus,, 
the intermolecular anisotropic potential for B^-X can be 


written as 
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[v A (r,8'3 H2 _ Ar = b l (r) P 2 <cos 6 X ) 

[v A (r,9- ! )] H = b 1 (r) [P 2 (cosa|_) + p 2 ^ cos9 2^ 


(3.72) 


? 2 
H2-N2O 


+ b 2 (r) j^a q Y 2 q^'l^ Y 2o} £l '2 ) 


(3.73) 


(where the terms have the meaning as given in eq. (3.28) and 

eq. (3.29)) 

where 

b 1 (r)=-^ 9 - ^ , x = r/a (3.74) 


12 " T 

x x 


b b 47rQ H Q X 

b 2 (r) = ~Y 2 ^ where b 2 = — — , X = N 2 or 

X X 3 - ^ 

x = r/a 


(3.75) 


As described earlier, H 2 ~Ar collisions cannot produce 
transitions in the J states of the ortho~H 2 molecule since 
Ar is a spherically symmetric atom. H 2 ~1T 2 and H 2 *-N 2 0 
collisions also cannot produce any J transitions due to the 
mismatch in the energy differences betxreen the J states of 
N 2 , N 2 0 and H 2< Evaluating Q 2 (J, J) etc. for ortho~H 2 and X 
interactions and proceeding in the same way as done for H 2 
for the case of infrequent transitions between J states, the 
relaxation rate, (p /Hj) ~ , due to ortho-H 2 and X collisions 
only, are obtained as 

p H 2 -X 1.02 x 10 13 <D( J)> 

( #T") = T J ■£ CraT 

1 p l 4,-X (2 *Pt*H 2 -X )2( V 0 A ) 2 


(3.76) 
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for X = Ug* NO or m °l ecu l e having dominant quadrupole 

moment and 


(mM 

X 1 


H 2 -z 


1.02 x 10 15 <D( J)> 

P 1 4 2 -X ( 2lc P^H 2 -X^ 2 ^ k o^ 2 ” 


(3.77) 


for I = Ar or any spherically symmetric atom. In eqs. (3.76) 
and (3.77)» D(J) has the same meaning as given earlier in 
eq. (3.67), a^ and _ x are the parameters (see eq. (3.49)) 
corresponding to H 2 -X pair and Cj- is a quantity defined in 
eq. (3.48) in which average is taken over the rotational 

9 

states of X. Values of C^- are obtained as 5/2 for N 2 and K^O • 
P is the density in amagats and = 2.69 x 10 ^ molecule/cc 
is the number density corresponding to one amagat. 


I II . 3 Experimental Results and Interpretation 
III. 3.1 H 2 -Ar “xx'k 1 ' 11 * 63 

The proton spin- lattice relaxation times T^ were measured 
in H 2 ~Ar mixtures as a function of density and composition 
in the temperature range of 300-600K. Typical plots of T^ 
versus p at different temperatures are shown in Pigs. 10 
and 11. 

The values of T^/p for different compositions and temp- 
eratures are shown in Pig. 12. The data did not show any non- 
linearity in T^/P versus composition and so (T^/p) 0- ^ 1 * 


values 



in milliseconds 



Fig. 11 T| vs density at 600 K in b^-Ar mixtures- 




300 K 
400 K 
500 K 



p6DUUD / SpUOD0SI ] | ltd U! 8/1 1 


Fig. 12 Dependence of T|/p on Argon concentration at different temperatures 
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for ortho-H^ and lr collisions alone were obtained by using 

a linear extrapolation of the plot to 100$ Ar. The solid 

line represents linear fits of the data using method of 

least squares. The error bars show the standard deviations. 

However, a non-linear behaviour for T / p versus composition 

20 

is predicted by uhe Bloom-Oppenheim two level theory 
which takes into account the transitions between different 
J states and assumes that only the rotational states J=1 
and J=3 for ortho~H 2 and the states J=0,2 and 4 for para-H 2 
are populated. This theory, however, cannot be used satis- 
factorily at temperatures above 450K since the fractional 
populations in the states higher than the first two for 
ortho-H 2 become appreciable. At present no theory which 
takes more than two levels into consideration including 
inter J state transitions is available. Hence to take into 
account any nonlinearity effects, if present at all, the 
T-^/P data were fitted to an empirical relation of the form, 

T±/P = A Q + A-^ + A 2 X x (3.78) 

where X^. denotes the fractional concentration of the compo- 
nent X in H 0 -X mixture, as has already been done for H 0 -He 

21 o 7 2 

mixtures by Lalita and Bloom . ( T^/ p) is obtained by 

putting X x = 1 in eq, (3.78). 

l p O A t* 

Poster et al. have obtained (T^/p) ~ values in the 
temperature range of 160-270K using a linear variation of ■ 
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T^/ p versus composition* These values as well as the values 

O ™ At** 

of (T^/P) obtained, by the linear and empirical fits above 

room temperature are shown in Fig. 13 as a function of temp- 
erature. From the data it can be seen that the difference in 
the values of (T-j/p)^ ^ obtained from linear fit and 
empirical fit is maximum at room temperature and negligible 
at 600K. Hence it may be inferred that the deviation from 
the linear dependence of 1^/p on composition is most prominent 
around room temperature. The values of (h o ) ( ^“'^ r are obtained 
from the experimental values of (^/p) 0 '"^ using eq. (3.77), 
and are shown in Fig. 14 and given in table 3.2 as a, function 
of temperature. These values are referred to as ( k 0 )g^pti 
and the corresponding (T 1 /p) 0-Ar values as ( ^p/p) gxpt] ' 


Interpretation 


The inter mole cul a r interaction potential between H ^ 
and Ar can be written in the following form (see eqs. (3.68), 
(3*69), (3.72) and (3.74)). 

V = ( r ) + v A (r,e • ) 


= 4e (• 


x 


1 

12 


3/ £L 

^g-) + (“Y5 ~ ~§) ^2 ( cos 6- 1 ) > x 

X XX 


a 


(3.79) 


where orientations are with respect to f (i.e. in molecule 
fixed frame whose Z axis is along r) . One can then write 
using eqs. (3.61), (3.77) and (3.51) 
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(6/LjJ 

jv-o 


}D6DUJD / 30SUU 
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(k o } 


0-£r 


^ ^ 2 ^H 2 -Ar 


= [ a^ 2 I ( 2,12) + a 2 a l(2,6)-2a 1 a 2 l(2 J 6,12}] (3.80) 

where 1(2,12) etc. are certain integrals defined, in eqs. (3.52) 
and (3-53). These integrals were evaluated on IBM 1401 and 
7044 computers at the I.I.T. Kanpur Computer Centre using 
Gaussian Quadrature Technique. Values of these integrals for 
H 2 ~Ar are listed in table 3.3* 

The constants and a 2 , which are the strengths of 
the repulsive and attractive terms in the in ter molecular 
anisotropic potential, were calculated from the above relation 
using the method of least squares. The values of a^ and a 2 
are oh-train©* tko values of ( corr CSp01lding 

to two different sets of (iq/^exptl o^sined from linear 
and empirical fits and are given in table 3.4. The solid 
lines in Pig. 14 represent the calculated values of ( k Q ) ^ ’"'^ r 
using the eq. (3.80) and the values of a^ and a 2 obt ain ed 
for linear fit. 


The relative anisotropy, a 2 /4e, in the attractive term 
of the intermole cular potential was obtained as 0.137 for the 
linear fit and as 0.135 for the empirical fit. it may be 
seen from the eq. (3.69) that a 2 /4e = ^-((k,, - a ± )/a) H . Using 
the polarizability data its value is obtained as 0.128 
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Table 3.2. Values of (i^/p) 0- ^ and as obtained by 

fitting the T^/o vs. composition data by a linear 
relation and by the empirical relation (eq. (3.78)) 


Temp. 

°Z 

T-jV P versus composition 

fitted b 3 r a linear 
relation 

1 ■ — 

| Tp/p versus composition 
j fitted by the empirical 
i relation 

(fp/P) 0- -^ 

in msec/ i 
amagat 

J 

(k )°" lr 
v o'exptl 

(l0 24 sec -2 ) 

(n/p) 0 -" 

in msec/ 
amagat 

(k )°" Ar 
^ o ' exptl 

(10 24 sec -2 ) 

164 

i 

0.082 

1.14 



185 

0.079 

1. 27 



204 

0.076 

1.39 

j 


225 

0.071 

1.53 



246 

0.065 

1.65 I 

j 



267 

0.061 

1.79 



300 

0.060 

2.21 

0.055 

1.95 

400 

0.052 

3.63 

0.048 

3.33 

500 

0.056 

6.05 

0.056 

6.02 

600 

0.060 

10.10 

. 0.059 

9.92 

i 
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Table 3.3. Values of 1(2,12), 1(2,6) and 1(2,6,12) for dilute 
H 2 ~Ar mixture including the first quantum correc- 
tion for a Lennard- Jones isotropic potential 



= T5.118°E, (a) = 

( M-) H Jr = 1-919 amu] 

3.169JP , 

(3e 

1(2,12) 

1(2,6) 

1(2,6,12) 

0.45516 

0.03427 

0.05823 

0.04H3 

0.41723 

0.03639 

0.05844 

0.04240 

0.37551 

0.03934 

0.05896 

0.04421 

0.34137 

0.04244 

0.05965 

0.04612 

0.31292 

0. 04566 

0.06045 

0.04809 

0.28340 

0.04979 

0.06156 

0.05064 

0.25720 

0.05458 

0.06287 

0.05350 

0.18590 

0.07602 

0.06867 

0.06575 

0.15170 

0.09537 

0.07348 

0.07602 

0.12520 

0.11941 

0.07891 

0.08799 
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Table 3.4. Values of the strengths ai and a 2 of the repulsive 
and attractive terms of the anisotropic part of 
the potential (eq. (3.79)) when fp/p versus composi- 
tion is fitted by linear as well as by empirical 
relation (eq. (3*78)) for H 2 -Ar. 


When ( T-j^/p) 0 Ar 

a 1 x 10 15 

a 2 x 10 15 

a 2 /4e 

is obtained by 

( ergs) 

(ergs) 


linear fit 

8.19 + 0.01 

5.68 + 0.01 

0.137 

Empirical fit 

8.11 + 0.01 

5.59 + 0.01 

0.135 
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In Pig. 14 the point corresponding to the (k )®~& r 

o'exptl 

value at 600K lies slightly higher than the calculated value 
using the a^ and obtained from the least square fit. In 
order to get a better fit of the entire data s anisotropic 
potentials as given in the eq.(3.79) but having repulsive 


—12 

terms with radial dependence other than r were tried 

‘""H. * 

(namely, r where n' =14, 15, 16, 17, 18 and 19. Higher 

values of n’ were chosen since the temperature dependence of 

the integrals l(p) becomes more steeper as n 1 increases). 

It was found that the model potential containing the repulsive 

“12 

term varying as r represented the data best. (These dif- 
ferent fits are not shown in Pig. 14 for clarity. ) y^g 
large difference in (k Q ) 0 “ Ar values at 600k might be due to 
the fact that in obtaining the values of (1^/ P )g~p^ for 100$ 
Ar in H 2 ~Ar mixtures a linear (or an empirical relation 
given by eq. (3.78)) has been used whereas a theory taking 
into account all the populated J states and the inter J 
transitions in a I^-Ar mixture might yield some other depen- 
dence of Tp/p- on the composition which should be used to 
oMain (V p >e^W- 


The value of a 2/4s obtained using Bloom-Oppenheim 
theory and the model potential given in eq. (3.79) i s close 
to the value obtained by polarizability data. Prom the agree- 
ment of the two values of the relative anisotropy, a^/4-z, in 
the attractive term of anisotropic intermole cular potential, 
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obtained from the px-esent analysis over a vri.de temperature 
range and the discussion presented above it may be concluded 
that a model potential written as 

r 8.11xl0~ 13 5. 

V 1 o ' 


V 


4e [■ 


1 

~12 


x 


X 


— — ^2 ( cos8 l) 


or, 

Y = 41.46x1c" 15 erg [(|) 12 -(J) 6 ] + 5 . 5SxlO~ l5 erg 


x [1.45(|) 12 -(|) 6 ] P 2 ( cos .9 -J_) (3.81) 

where a = 3.169 Jl for H 2 *-Ar pair represents the MER data well 
and also that the polarizability data can be used to estimate 

■H 

the strength of the attractive anisotropic inter molecular 
potentials. 

IP 

Poster et al. analysed their T^ / P data in H 2 ~Ar and 
other H 2 ~inert gas mixtures at low temperatures using Quantum 
Scattering theory and assuming an intermole cular potential of 
the form as given in eq. (3.79). They obtained values of the 
relative anisotropy, a 2 /4e, in the anisotropic attractive 
potential which are twice of that obtained by polarizability 
data whereby they suspected that perhaps the radial dependence 
of the repulsive part of the anisotropic potential might be 
other than r 12 . But Hie hi et al. 11 have fitted the H^-He 
data in low temperature range using an exp(-6) potential and 
obtained a value of a 2 /4e which is again twice of that 
obtained by polarizability data. So it may be concluded that 
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the Quantum Scattering theory does not give a value for 
a 2 /4s consistent -with polarizability data. 

I II . 3. 2 H 2 -N 2 Mixtures 

1^ in H 2 “-IT 2 mixtures were measured a.c a function of 
density and composition in the temperature range of 300-6002. 
The plots of T-^ versus p and T-, / p versus composition are 
shorn in Pigs. 15, 16 and 17. The solid lines represent 
linear fits. The error bars denote the standard deviations. 
Where these are not shown, the standard deviations are less 
than 2/£. 

In this case also non-linearity of T^/ d versus compo- 
sition was not observed and thus values of (T,/p)°“ N 2 were 
obtained by using linear extrapolation. But the two level 
theory and a recent work by Kalechstein and Armstrong^ 
predict a non-linearity. Hence ( T^/p ) values were also 

found out by fitting the data to the empirical relation as 
given by eq. (3.78) in a manner described for H 2 ~Ar mixtures. 
It was observed that the values of ( obtained from 
the empirical relation are very close to those obtained by 
a linear fit within the experimental error. These values 
of (T^/p)^~®” 2 together with those obtained by Williams^ 
in the temperature range 150-2952, which they have obtained 
by assuming a linear variation of T-^/p versus composition 
are shown in Pig. 18. Williams' data were also fitted to the 
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empirical relation. The values of (^/p) 0 "^ 2 are shown in 
Fig. 18. All these values of ( T^/p) < “ > "”^ 2 were used to calcu- 
late (k Q + Cjj k 1 ) C " J "2 using ea. (3.76). The values of 

and, of (k o +C^ k- L )^“'^2 are also given in table 3.5. 

2 on 

Fig . 19 shows the values of (k o +C N 1^) "’'2 obtained from 
the values of (T^/p) 0 ^" 2 evaluated by a linear fit of T,/p- 
versus composition data. The error bars represent standard 
deviations . 


Interpretation 

The intermole.cular potential between Hg and IT^ Can be 
written as (see eqs. (3.68), (3.70), (5.73) and (3.75)) 


V 


= + (~y^ §)[P 2 (cos 0 ^) -i- PgCcose’g)] 

x x xr~ x ~~ 


+ ( ^2 " ^ 5 |_2 Vtq^P r 2q (a 2 } 

with X' = r/a. 


l3 2> 2 


(3.82) 


For this form of potential one can write, using eqs. (3.76), 
(3*61), (3,62)- and (3.51) 


0-IL 


(k o +C N„V * = I ^ 2 » 12 )' + a 2 l(2,6):-2a L a 2 l(2,6,12)] 


+ CL ( — ™)[bj I(4,12)+b!l(4»5)-2b b 1 ( 4 , 5 , 12 )] (3.83) 

2 5 0 7t.il 2 x 


where CL. = 5/2 and I (4,12) etc are certain integrals an 
iM 2 

defined in eqs. (3*52) and (3.53). The values, of the integrals 
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0-N o 0-N o 

Table 3.5= Values of (i^/p-) 2 and (k^Cfr h^) 2 as 

obtained by fitting the f-^/p versus composition 
data by a linear relation and by the empirical 
relation (eq. (3.78))., 


Temp. 

( °K) 

T-j^/p versus composition 
fitted by a linear relation 

Tj/o versus composition 
fitted by the empirical 
relation 

( Vp ) 0-1 * 2 

(msec/amagat) 

< Wi >°' K2 

(10 2 5 see -2 ) 

( V d °~ lf 2 

( msec/arnagat) 

(VQlah ) 0 '® 2 

(10^5 sec" 2 ) 

150 

0.860 

H 

H 

• 

H 

0.822 

1.07 

170 

0.793 

1.15 

0.785 

1.14 

200 

0.687 

1.20 

0.680 

1.19 

295 

1 0.500 

1.62 

0.520 

1.68 

300 

0.510 

1.66 

0.540 

1.76 

400 

0.330 

1.95 

0.350 

2.06 

500 

0.215 | 

1 2.05 

0.212 

2.00 

600 

0.170 

2.47 

0.163 

2.39 
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Table 3.6. Values of 1(2,12), 1(2,6) and 1(2,6,12) for dilute 
H 2 -N 2 mixture including the first quantum correc- 
tion for Lennard-- Jones isotropic potential. 



[UA) H2 . 

~N 2 ~ 59 *44 ^ a ^H 2 -lT 2 - ^* 

314A°, 


^h 2 -u 2 

= 1.881 amu]^® 


|3£ 

1(2,12) 

1(2,6) 

1(2,6,12) 

0,39624 

0.03780 

0.05867 

0.04326 

0.34962 

x o 

H 

^r 

O 

• 

O 

0.05946 

0.04562 

0.29717 

0.04779 

0.06101 

0.04940 

0.20148 

0.06975 

0.06702 

0.06227 

0.19812 

0.07101 

0.06735 

0.06297 

0.14712 

0.09879 

0.07429 

0.07777 

0.12007 

0.12554 

0.08021 

0.09094 

0.09906 

0.15899 

0.08678 

0.10631 
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Table 3.7. 

Values of 1(4,12), 1(4,5) and 1(4,5,12) for 
dilute H 2 -N 2 mixture including the first quantum 
correction for iennard- Jones isotropic potential 

Pe 

1(4,12) 

1(4,5) 

1(4,5,12) 

0.39624 

0.02757 

0.05015 

0.03170 

0.34962 

0.03031 

0.05018 

0.03314 

0.29717 

0.03479 

0. 05066 

0.03550 

0.20148 

0.05059 

0.05330 

0.04357 

0.19812 

0.05150 

0.05346 

0.04399 

0.14712 

0.07155 

0.05686 

0.05316 

0.12007 

0.09087 

0.05990 

0.06143 

0.09906 

0.11463 

0.06322 

0.07061 
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are given in Tables 3.6 and 3.7. 

Prom the eq. (3.76) it can be seen that the experimental 

observation of the proton spin-lattice relaxation tines T‘ in 

Hp-KL mixtures gives information on (k + C. T v ) 0-W 2 g^d not 
^ ^ o 1' 

on k Q and k^ separately. Therefore to extract information 

about the attractive and repulsive parts of the anisotropic 

in ter molecular potential (eq. (3.82)) the (k +0^ k 1 )°“ IT2 

values were fitted to eq. (3.83) using the method of least 

squares and considering the cases where k is assumed to be 

much smaller than k^ and hence neglected as was done for 

2 21 

H 2~ C0 2 mix ' Jc ures by previous workers and when the effect of 
k Q is also taken into account. The values of b^, b^ and 
a l s a 2 y ^i» ^2 * or ^ ese f w0 cases were calculated. Neither 
of these two models gave satisfactory fit of the experimental 
data. Hence slight variations of these models were tried in 
which b 2 was taken to be equal to the strength of the quadru- 
pole-quadrupole interaction between the molecules and 
The repulsive term in k Q was excluded and the constants a p , 
b-^ were evaluated. This also did not give a good fit of the 
experimental data nor it gave a value of a^ consistent with 
the polarizability data. Another model was tried in which 
the repulsion in k-^ was excluded and a^ and a^ were evaluated 
using the method ox least squares. It was found that this 
model represented the data best as can be seen from the 
Pig. 19. a-j_ and a 2 for this model are obtained as 
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^ -j r~ 

( 6 . 98+0. 01) xlO ergs and (4.25+0.01)xl0 - ergs respecti- 
vely. The relative anisotropy, a 2 /4e, in the attractive r^ 
term of the anisotropic part of the potential was obtained 
as 0.131 which is close to the value, 0.128, obtained from 
polarizability data. 


Thus it may be concluded that a potential written as 
V = 4s (-^2 ~ 4^25x10 ) [ p 2 ( cos6 / ) + P 2 ( c °s q^) ] 

IK X X X 

4tcQh 2 Qu 2 ° 


5 3 ’A V a l Y 2q ( ^V ^q^ 2 ) 


a x 


q=-2 


or, 


Y = 32.81x10 15 erg [ (~) 12 -(§ ) 6 ] + 4.25xlO~ 15 erg[l. 64(f) 12 -(f ) 6 ] 

x [P 2 ( cos 6^) + P 2 ( cos QbJ ] 


27.73xl0“ 15 erg (f)^ * J ^ a q T 2q (ft{) 2 ) (3.84) 


a\ 5 


where the values e = 8.203 x 10 ergs, a = 3.314 J? for H 2 -H 2 

"*■26 

pair for lennard- Jones potential and Q T - = 0.63 x 10 e.s.u. 

? -26 9^52 

cm , and Q. T = 1.4 x 10 e.s.u. cur are used, 

2 

represents the IT.M.R. data well and also that the polarizability 
data can be used to get an estimate of the strength of the 
attractive r 6 term of the anisotropic intermolecular poten- 


tial. 



88 


III .3. 3 H 2 -N 2 0 matures 

Proton spin-lattice relaxation times T^ were measured 
in H 2 -N 2 0 mixtures as a function of temperature, density and 
composition in the temperature range of 200-400K. The measure- 
ments were done up to 400K only since some chemical reaction 
was observed to occur in the mixture inside the Be-Cu pressure 
vessel at higher temperatures. Typical plots of the values of 
as a function of p- and of T,/p>aS a function of composition 
are shown in Pigs. 20, 21 and 22. 

In this case also T^/ p was found to vary linearly with 
composition. The analysis was done exactly in the same way as 
was done for mixtures. It was observed here also that 

values of ( T^ / P) ^ ^2® obtained by empirical, and linear fits 
are close to each other (see table 3.8) and thus one can consi- 
der T^/ P to be a linear function of composition. Values of 
( T^/ p) stained using a .linear fit and the empirical fit 

are shown in Pig. 23 along with the data obtained by Johnson and 
Waugh"’ at room temperature. The values of (k Q +C„ 
were calculated using these values of (T^/bO^ - ^ 2 ^ and eq. (3.76) 
and are given in table 3.8. Values of (k 0 +C^ q^) O - ^ 2 ^ 
corresponding to the (^/p) 0-1120 obtained by linear fit are 
shown in Pig. 24. 

Interpretation 

The intermole cular potential for H 2 and molecules 



65-4 %N 2 0 


>^^19.5 %N 2 0 


)’ ”10 20 30 40 50 “60 

Density in amagats 

Fig. 2 0 Plots of Tj vs p in H 2 -N 2 0 mixtures 


65-4 %N 2 0 


19.5 % N?0 


Density in amagats 

Fig. 21 Plots of Tt vs p in H 2 -N 2 0 mixtures 
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Fig. 23 Values of (T|/^)° as a function of 
temperature . 
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can be written,, neglecting, the small dipole moment on , 

in the form of eq. 3. 82, . whence (k 0 +C N qIc,) 0- ^ 0 can be 

written in the form of eq. (3*83). The values of the 

constants a-^, a.^ b-^ and b£ were calculated using the method 

of least squares. (The values of integrals 1(2,12) etc. 

occurring in eq. (3.83) are given in Tables 3.9 and 3.10.) 

Different model potentials were tried to fit the data best, 

as was done for TL--1L, for the cases when k ■ << C. x n k, and 

2 2 o IT^O 1 

when k Q -w q^. It was found that a potential of the form 
of eq. (3.82) with 


= 0 


> 


a. 


0 


r e - 


k 0 << 


represents the data best as can be seen from the Fig. 24 
(fits with other forms are not shown in this figure for 
clarity), b-^ and obtained in this way are (10. 33+0.0l)xl0" 
ergs and ( 4.93+0.09)xl0~ k ^ ergs respectively. The model poten- 
tial (eq. (3.82)) 'then can be expressed as 


V = 4e [-ij - 

X X 

. ,10.33xl0 -14 4.93X10 -14 -, y „ ,r L ' , ,,js n ' , 

+ (--^5 -5 ) 

erg 

(3.85) 

with x = r/a and e = 11.54x10 ^ ergs for ^2"~^2 G 
quadrupole moment 


of N 9 0 was obtained using the relation 

■nU c. 


■14 
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Table 5.9. 


Values of 1(2, 12), 1(2,6) and 1(2,6,12) for 
dilute H 2 “U 2 0 mixture including tlie first 
quantum correction for Lennar cl- Jones isotropic 
potential. 


[ (e/lO^^o = 85. 6 24° K, (a) Il2i . ll2 o 

^Ha-NaO = 1.928 amu. ] 48 . 


5.759a, 


pe 

1(2,12) 

1(2,6) 

1(2,6,12) 

0.41812 

0.05627 

0105855 

0.04242 

0.57166 

0.05971 

0.05907 

0.04445 

0.55450 

0..04522 

0.0 '5988 

0.04662 

0.27875 

0.05080 

0.06196 

0.05150 

0.24050 

0.05864 

0.06422 

0.05602 

0.20648 

0.06807 

0. 06662 

0.06156 
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Table 3.10. 

Values of 1 ( 4 , 12 ), 1(4,5) and 
dilute H 2 -N 2 O mixture includin 
quantum correction for Lennard 
potential. 

1(4,5,12) for 

2 the first 
-Jones isotropic 

Pe 

1(4,12) 

1(4,5) 

1(4,5,12) 

0.41812 

0.02658 

0.05030 

0.03124 

0.37166 

0.02890 

0.05013 

0.03239 

0.33450 

0.03139 

0.05026 

0.03371 

0.27875 

0.03685 

0.05098 

0.03660 

0.24030 

0.04235 

0.05 188 

0.03944 

0.20648 

0.04934 

0.05309 

0.04295 
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4 * Qh 2 %?0 


b 2 " a 5 

a H 2 -lT 2 0 


CL. = 0.6 3 x 10" esu cm 2 

lip 


a H 2 -N 2 0 


3.155 


as (4.65 + 0.09) x 10 26 esu cm 2 . I his value is higher than 
the recommended value which is 5.0 x 10~ 2 ^ esu cm 2 ^ 2 


though it is in agreement with the value obtained earlier 
by NMR method (Q Q = 4.25 x 10~ 2 ^ esu cm 2 )^ using the data 
at room temperature. It is quite likely that the form of the 
repulsive part of the anisotropic potential used here needs 


to be improved. 



CHAPTER IV* 

PROTON SPIN-LATTICE RELAXATION IN CII^Cl GAS 


IV.l General 

CH^Cl is a symmetric top molecule having C^ v symmetry. 
Its rotational states can be described by the quantum numbers 
J, K, and M which correspond to the total angular momentum 
the projection of J on the symmetry axis of the molecule and 
the projection of J on a space-fixed Z-axis respectively. 

If denotes the rotational Hamiltonian for a symmetric 

top molecule, one can write 

%j |JKM> = fi m. jE | JKM> 

where 

= I T~ J ( J+ D + (y~ ~ r~) k2 (4-1) 

o <■ ± a o 

= hBJ (J+l) + h (A-B) X 2 (4.2) 

where 1 ^ and I Q are the moments of inertia of the molecule 
parallel and perpendicular to the symmetry axis and their 

55 A paper entitled ’Proton spin-lattice Relaxation in Low 
Density CH 3 CI Gas* based on the material presented in this 
chapter was published in Chemical Physics Letters, Vol. 58, 
page 575-378 (1978). 
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values for CH^Cl are^ 5.49 x 10”^ gm cm 2 and 57 x 10 ^ 
gm cm 2 respectively. In 'writing the eq. (4.1) it has been 
assumed that the molecule is in ground electronic and 
vibrational states. 

IV. 2 Review of tho Spin-relaxation Th eory 

The dominant intramolecular interactions that contri- 
bute to the spin-lattice relaxation in dilute gases are the 
intramolecular dipolar and spin -rotational interactions. 
Assuming that the contribution from the intramolecular 
dipolar interaction is negligible for CH-Cl as was founa for 
many symmetric top molecules BP^» CHP-, CII^P, HU by 

bong and Bloom’ ^ and Armstrong and Courtney , the proton 
spin-lattice relaxation rate can be written as (see eq.(l.22)), 

= R C 

The theory of nuclear spin relaxation due to spin -rotation 

interaction was first discussed by Hub bar d^ and Bli char ski ^ 

43 

and was extended by Bloom et al. for the case of dilute 
polyatomic gases. They have obtained expressions for T^ in 
terms of the Pourier transforms of the correlation functions 
of the spin-rotation interaction. A brief review is presented 


below. 
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The Hamiltonian describing the spin- rotation interac- 
tion can be written, in general, as 

$. s . r> = -hf-C-J (4.3) 

where I is the nuclear spin angular momentum, J is the 
molecular rotational angular momentum and 6 is the spin- 
rotation tensor. In polyatomic molecules, the spin-rotation 
interaction is in general characterized by six independent 
parameters, viz., the three principal values of the spin- 
rotation tensor and the three Euler angles required to 
specify the orientation of the principal axis coordinate 
system. Eor symmetric top molecules of symmetry (such 

as molecules of the form YX^ or CX^Y) the spin-rotation 
tensor for the I nucleus which is on the body sy mm etry axis 
is of the form 



where Cj_ and C„ are proportional to the magnetic fields per 
unit angular momentum produced at Y site by rotations of the 
molecule about axes perpendicular and parallel to the symmetry 
axis respectively. Thus, only two independent parameters, 
namely, and Gj t are needed to characterize the spin- 
rotation tensor. It is customary to represent the spin- 
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rotation tensor in terms of quantities C o and Ch defined 

' Co 

as 


C 

a 



+ 2Cj_) , 



(4.5) 


C is the spin-rotation coupling constant for the isotropic 
a. 

part of the interaction whereas is for tlie anisotropic 
part of the interaction. The Hamiltonian for the spin -"rota- 
tion interaction for the Y nuclei may he 'written as ^ 


ii- 1 ^ 


s.r . a 


0. . f C" 1 )* liJl k -C d do J o- I -li- 1)kl i J -0 


k=-l 


(4.6) 


where 


I„ = I. 


L + 


+ ( X + i I ) 

If 2 v x - y' 


and 


J_ = J. 


J + = 


+ fj ( J z t 1 V 


The primes in eq. (4.6) denote that the components are taken 
in molecule-fixed coordinate frame with its Z axis along the 
molecular symmetry axis. For one of the off-axis X nuclei, C 
can he written as^'*^ 


C 


0 

0 sx" 0 y y 

0 


c xx C yy 
0 
0 


0 

0 

0 


H- 
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^ C zx- C y P 

0 

0 


0 

“X C xx” C yy^ 


0 

0 

0 


+ 


^ C xx +C yy^ 
o 


0 


0 


^■ (c xx +c yy ) 


o 


o 

o 

C, 


ZZ 


(4.7) 


where G C arid C are the diagonal components of the 
XX yy zz 

spin-rotation tensor for one of the X nuclei. The first two 
terms have matrix elements connecting the states for which 
AK = +1, +2 and the last term for AK = 0. For dilute gases, 
it can he shown^’^ that the matrix elements between dif- 
ferent K states can be neglected. C, then, contains only the 
last term in eq. (4.7) and only the parameters C A and G\j 
defined as, 


and 


C 


1 


= 2- 


^ ^XX + ’X-cr) 


yy 


C.i 


zz 


(4.8) 


are required to describe the spin-rotation tensor. The spin- 
rotation tensor for X nuclei also, thus, reduces to the form 
of eq. (4.4). The Hamiltonian for the spin-rotation inter- 
action for the X nuclei, therefore, can be written in the 
form of eq. (4.6) where C and C-, would be given by eqs . 

(4.8) and (4.5). 
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In space fixed frame the eq . (4.6) becomes 




s 


r = S' Z m *1 

• m 


m 


(4.9) 


where 


K 


lm = -(f) J 0 d I (-l)^0(112 !m(1 )D o 2 > _ (iw . ti) J j (4.10 


P- 


in which D are the rotational matrices (Rose^) . 

Using the general theory' 1 ', R^ is written as, 

R C ~ 4%2 J lm 

where J lm ( w q ) is the Fourier transform of U^m^ 


(4.11) 


CO 


( 4 . 12 ) 


m 


(4.13) 


- -lfcLT 

y%) = J hm {x) e dT 

"“co 

is the correlation function of 

<W T > = K e <E i m (0) 

where the symbol < > represents an ensenible average. For a 
freely rotating molecule [ U lm ( t) ] f ree comes out to be 
in ae pendent of m. Hardy ^ has Shown that -under quite general 
conditions, one can write 


G lm (T) , = U 10 (t) exp (im 60 jt) (4.14) 

where ta-j is the Larmor frequency for the molecule. Thus, in 
order to evaluate G- lm (m) only & 10 (t) need to be evaluated 
which, with the help of eq. (4.13), can be written as. 
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g io ( t) 


rJc 2 < J o(°) «L(r)> 

e L a o v ' o ' ' 


2(f) *C a C d l '(-1)^C(112jO|_i) <<_ U (0)J (0)J o (t)> 


P- 


+ | C 2 I _ (-l)' a+ ^’c(112 5 0|a) 0(112,0^') 


P>P- 


* i 3 l~4 0) I>^_ u ,(t)> ] (4.15) 


P 


P 


P 


Due to collisions between molecules the time dependence of 
G’IqC'O is affected. for polyatomic molecules it was assumed 
by Bloom, Bridges and Hardy ^ that the effect of collisions 
could be adequately described by writing 


^im^-free (4.16) 

where g 2 ( 0 ) =1 and g 2 (i;) is a function which decreases 
monotonically to zero as % is increased and is called the 
reduced correlation function. So, to take into account the 
effect of collisions on and hence on G^q all the correla- 

tion functions occurring in the right hand side of eq. (4.15) 
are expressed as products of corresponding free molecular 
correlation functions with monotonically decreasing functions 
of time g 1 (r), g ? (r) and g^ 2 (r) respectively, i.e. 


<J o (0)J q (t)> 

y°) 


= t <J O (0)J o (T)>] free g l( x) 
j 0 (x)> = [<1^(0) jyo) J 0 ( 


T ) > -^free 


(4.17) 

s 12 (t ) 


(4.18) 
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<<-4°) V 0) V (t) J £-h' ( t)> ' 


= r< D o,-l°> V 0) j u' ( t) 


free *=12 ^ 


0>-pr pr - p, 1 

(4.19) 

where g-j^r) , S*2^ T ^ afth. ^12^ are "^he 1 reduced correlation 
functions' with g-, (0) = g^( 0) = g^ 2 (0) = After evaluating 

the free-molecular correlation functions it is seen that 
they contain terms that are time independent as well as 
terms that oscillate at frequencies corresponding to rota- 
tional transitions AJ = +1. For dilute gases it can he shown 
that the oscillating terms do not contribute appreciably to 
the spectral density at the larmor frequency and, hence, can 
be neglected, 'i'hen with the help of eqs. ( 4 . 17), (4.19) and 

(4.15) one get s 


g 10 (t) = i C 2 


OQ 


3 a 


r 1 ( 2 j+i) j( j+i) p 


J=0 E=-J 


JEM 


xfspv - ftr [2^^ 


(x) 




r 3K 2 -j( J+l)-, 2 

t—JTJ+iT ] 


g^ 2 (T)} 


(4.20) 


fhe spectral density J 10 (w) is obtained using eqs. ( 4 . 15) and 
(4.12) as 

J 10 M = P(») 5 C 2 a <J( J+l)> ^(m) 


( 4 . 21 ) 
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where 

and 

P (.“.) 


3 l(a>) = J. e 


-10 T 


dr ; C0=.iO 0 -tOj 


9 C, 

i-4 (^-o 

^ a 


< 3 K 2 - j( j+iy - ^12 
'7jTj+iT> ■ 


( tl) 


3^TuT)~ 


, i , c d,2 <jrmr - “ 2 + j( J+D* 

9 l cj ' ' <J Cj+iT>' 

ct 




(4.22) 


and anc ^ Dqo^) are spectral densities of g^^T) 

and g^(x) respectively. Using eqs. (4.21), (4.22) and (4.11) 


and assuming that . ai Q >> ■ -a>j , is obtained as 


H 0 = if- p(« b ) Cl <J(J+1)> }l (a 0 ) 


(4.23) 


If all the 'reduced correlation functions' occurring in eq. 
(4.20) are assumed to decay with the same time constant x-^, 
i.e. 


g(x) = exp (-x/xj) 
then = d 12 (o>) = ij_ 2 ( w ) 

defining quantity a as, 


(4.24) 


2^1 

-i _i_ 2 2 

-L 1 ♦ (a) % 


a 


■fu 


2I 0 kT 


(4.25) 


and evaluating the averages and using eq. (4.22), the eq. 
(4.23) reduces to^ 


R o = 


2m 4 


[Ini 0 ~2(I 


a '•"Ola 


10 3 01' ad 

T 1 


I m) c „ c ri + ( I 2~l~ 3*10 + 9 x 01^ C d^ 


“I i 2 2 

l+o) 


(4.26) 
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where 


where 


-01 


-10 


■ 2-1 


• 1-1 


2a 

3“ 

"41 ( J+l )> = 

2a 

<E 2 > = ; 

3 


2a 

< . E 4 

3 

jTj+i) 


1 -JL. 


(i=?-)Vl - S- + a) 

y - 1 


jr 4 - 

j( j+Ty 


> - 


= - (Mr^ )(1 - s) 


y 


y = 


S = 


(1 _ i4) 2 _ / P 7 

I o ~ ( «V 

2y ln for p > 0 


'i-y ; 

1 , _i 

- jy|' can | y | for p < 0. 
Eq. ( 4 . 26 ) can he written, in short, as 

- 2 TU 


R, 


-4jl_ P 2 

a 0 


eff 777“' 
1 + (w 0 


-x 


»j) M 


where 


(4.27) 


(4.28) 


c eff ~ (1 ” 3~'^ C 1 + f y 2 c a c d 

+ | [S(^f-) 2 (-l - + s ) _ i + | y 2] c d (4.29) 


The spin-lattice relaxation time due to spin-rotation 
interaction is then given by 
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1 _ An 2 n 2 

T ± ~ a " °eff 


1 + - w.j) ^ 


(4.30) 


IV. 2 Results a nd I nter lire tat ion 

i'he proton spin-lattice relaxation time was measured 
in CH^Cl in the -density range of 0. 2-6.5 amagats at room 
temperature ( 298K) . i'he dependence of on the density is 
shown in Rig. 25.' long relaxation times coupled with low 
signal to noise ratio made it difficult to take data below 
0.2 amagat . The details of sample purification and prepara- 
tion are given in Chapter II. 

i > 

Assuming that the spin-rotation interaction is the 
dominant mechanism of relaxation in CH^Cl the relaxation 
rate can be given by eq. (4.30). Prom this equation it is 
seen that the plot of T-^ versus P goes through a charac- 
teristic minimum at cot = 1 while for tdz 1 » 1, ^ 

Using the kinetic theory the correlation time may be taken 
as 

T 1 ~ p<0V> (4.31) 

where p is the number density of the gas and <ov>. is the 
velocity averaged collision cross-section for molecular 
reorientation. Since is inversely proportional to p, 
it is obvious that in the region cot 1 » 1, decreases as 
the frequency of collision ( arid hence p) increases but in 
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the region cOt-^ « 1, 1^'°= P (phenomenon called as pressure 

narrowing). Brom eqs. (4.30) and (4-31) and following 
45 

Armstrong et al. T-^ can he written as 


with 



A + B/ 


1 


a h 




'eft 


5nin 


, eo 0 » 


(4.32) 

(4.33) 


B 



1 - 

,2 “o '•“'min 

eff 


(4.34) 


where is 4he density at which T-^ minimum occurs. 

The experimental data between 0.2 to 6.5 amagats was 
fitted to eq. (4.32) using the method of least squares. The 
fitted curve is shown in Big. 25 "by solid line. The values 
of A and B obtained are given below along with their standard 
deviations 

A = 371.12+4.09 msec/amagat 

B = 27.55 + 4.03 msec.amagat. - 

The values of PJ^ and min were obtained as 0.273 amagat 
and 202.23 msecs, respectively. 

Using eq. (4.33) and the values of A and P min> the 
value of was obtained as 0.1 MHz neglecting in 

comparison to-w 0 . 
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The minimum value of T ^ occurs at (u 0 - = 1 


Since = P <crv> ' at constant temperature 


«tv> = 


u 0 “ 


min 


(4,35) 


Defining an effective cross-section, a e f-p> for molecular 

— _ ± 

reorientation as <tav>/ v, where v = (Qk'i'/jcp) 2 in which k 
is the Boltzmann constant and u, is the reduced mass of CH^-Cl, 
it was found that g Qff = 4.01 i 2 . Comparing tills value 
with the kinetic crosssection, g ^. n = % &^ T - where the 

A? 

value of the Lennar d- Jones parameter a CH C1 for methyl 
chloride is 4.15 2 . one get ^hin^eff *’ 14 , which indicates 
that it takes about 14 collisions on the average to randomize 
the rotational angular momentum. 



CHAPTER V 


ANALYSIS OF T 1 DATA IN HC1, *2 ' 1IH 3 ^ BP 3 


V.l Gene ral 

The general theory of spin— relafa ^io n expieoses oho 

relaxation rates in terms of the correlation inactions of the 

10 

intramolecular interactions. Using Bloom-Oppenheim theory 
the correlation functions for intramole cular intei actions can 
be written in terms of the correlation functions of inter- 
molecular interactions. Rigorous expressions for the relaxa- 
tion rates for simple systems like H 2 » nhere the number of 
rotational states populated at ordinary temperatures is one 

or two, were obtained by Bloom et al. 10 2X1 ^ " b3ae da ' ta 
were interpreted to obtain information about the anisotropic 
intermole cular potential. For systems of diatomic and 
polyatomic molecules where the number of populated rotational 
states is large, no such theory is available in literature. 
Hence, Rajan et all^ proposed a model according to which 
the correlation time corresponding to a s’tabe for the intra- 
molecular interactions is taken to be equal bo the lifetime 
of the molecule in that J state. The life bime is calculated 
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using the first order perturbation theory assuming that the 
collisions are weak. They applied this model to obtain the 
intermolecular potentials in spherical top molecules, CH^ , 
and SiP^ and in mixtures of CH^ with lie, he, Ar, and 
COg. Prom the analysis of T^ data, they obtained values of 
hyperpolarizability of CH^ and the octopole moments of CE^, 
CP^ and SiP^ which are in good agreement with the values that 
are obtained from other methods. Hence, one can conclude 
that this model worked well in these systems. 

The purpose of the present investigation is to test 
the validity of this model in some linear and symmetric top 
molecules. Two linear molecules having dipole and quadrupole 
moments, viz. HC1 and P^ respectively and two symmetric top 
molecules having dipole and quadrupole moments viz. , NH- 
and BP^ were chosen for this purpose. Por all these mole- 
cules, the spin-rotation interaction is the dominant 
mechanism of relaxation 

The Hamiltonian corresponding to the intramolecular 
spin-rotational interaction for the nuclei can be written, 
in general, as 1 

3o 

j'\.. = -hi • C • J (5.1) 

where I and J are the nuclear spin and molecular rotational 
angular momentum operators respectively and C represents the 
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spin 


in— rotation coupling constant which is a tensor for sym- 


metric top molecules and a scaler for linear molecules. 


Using the general theory, is given by 


(5.2) 


whe re 


~~ C 2 J(J+1) t 1 ' 


(5.3) 


for linear molecules, in the short correlation time limit. 
Here t, ^ is the correlation time for the spin-rotation 
interaction corresponding to the J state. Por a symmetric 
top molecule eq. (5.2) reduces to 


4JL_ c 2 t 
a u ef f T 1 


(5.4) 


where a and C 2 ^ are defined earlier in eqs. (4.25) and 
(4.29), respectively, in Chapter IV. 

17 T< 

According to Rajan’s model for a state J'M* 

can be written as 


2J'+1 


ll l . A( JM; J' M* ) 


J M=— J 


(5.5) 


M^M ; for J=J* , 


for linear molecules and 


( 2 J' +1) 2 M’ 


W W 

I • I 't 

=-J ! J M=- t 


MVM for J'=J 

for symmetric top molecules. 


I - A( JKMj J' K' M* )i (5.6) 

iac ( 

J 
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Here and A( JEM; J' K’ M' ) are the transition rates 

for the transitions J’M' JM and J'E'M' JEM for the linear 
and symmetric top molecules respectively and are evaluated 
in terms of the intermolecular interactions using first order 
perturbation theory (weak collision approximation) „ 

The intermolecular potential between two molecules can 
be written as 


V = v Q (r) + v A (r y e) (5.7) 

The isotropic part, V (r), is assumed to be a Lennar d-Jones 
potential 


v 0 (r> = 4e [(f ) 12 



e and a are the Lennard-Jones 
parameters , 

(5.8) 


and the anisotropic part, v A (r,9)„ c.an be expressed in terms 
of the multipole interactions between them.' 


The anisotropic attractive interaction potential 
between two molecules having axially symmetric charge 
distribution is ; in general, written as^^ a summation of terms 
of the form 




4* ^ (r) * l¥ , %ll2 m) ^(Sip 


= c 


W 

X r k+l 


(5.9) 


wi th 


^ < r > 
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where -2~ and. ft- are the orientations of the vectors 
shown in Fig. 26 related to a space-fixed frame (SFF). In 
this figure e^ and e ^ are elementary charge distributions 
on the molecules 1 and 2 respectively. Ihe molecules have 
been considered as having rigid non- over lapping charge 
distributions centered around origins 0-^ and 0^. The vector 
connecting 0-^ and 0^ denoted by r = (r,& -) where ^ = (qV.^') ? 
is along the in ter molecular axis. A = ^ = ^-q + ^ ^1 


and indicate the order of the electric moment z for dipole 
& 2 . = 1, the quadrupole 1 ^ = 2 etc. 

£ 


C, = 




( -1) r A-n ( 21 +1 ) i i 

! 2 - Taw l Ta^ryrTayTTT J 


1 

T 


( 5 . 10 ) 


The quantity Q denotes the scaler magnitude of the multi pole 

fl 

moment of order f or the axially symmetric charge distribu- 

tion. C m^m^m) is a Clebsch-Gordon coefficient as 
defined by Rose"^ and m are the spherical harmonics of 
order Z of the corresponding angles. 

The repulsive part of the anisotropic potential was 
-12 

assumed to vary as r having the same angle dependence 
as that of the attractive part. 


V. 2 Linear Molecules 

1 . Molecules having: dipole moment 

The anisotropic part of the inter molecular potential 
between two molecules having permanent dipole moments, such 
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as HC1, can toe written using eq. (5.9) , neglecting higher 
order interactions, as 


v A (r A e) 

= (ff) 2 to(ar) 

E.. C(112 J m 1 m m) Y-, ^ ( Q q) 

m l m 2 ' L 




(5.11) 

where 

to(a^) = 

x ^ 

a 2 

- -5 

X 

(5.12) 

with 

x = r/a and 

&2 = 4iTp, 2 /a 3 

(5.13) 


'the effect of the anisotropic potential is to produce 
transitions in molecule 1 (e.g. , HC1) from the state J',M’ to 
the state J,M while the molecule 2 (e.g., another IIC1 molecule) 
simultaneous undergoes a transition from J m , M m to J", M'\ 
Denoting the transition rate for this process toy 
A[ (• J 1 M' ) ( J"M"$ ) ] , one can write'*"*' 1 the transition 

rate for the transition J',M' ■+ ■ JM for the molecule 1 as 
(see eqs . (5.37) and (3.38)). 


A(JM 5 J'M' ) = £ . A[ ( JMj J' M 1 ) ( J"M n } ) ] (5.14) 


where 


= l. (ff) dfUjjjjMjn,) 


m l m 2 


x [C(U2 5 m 1 m 2 m)] 2 P JtI1M „, I ( JM I ^ ^3^ ) | J’M’) | 


| ( J"M n | Y lm (« 2 ) I ^"M"*) | 


X 


(5.15) 
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where ^ w jj'. ~ (Ej-Ej, ) + (Ej„ - Ej„, ) 


( 5 . 16 ) 


jm 




2J Tn + l’ ' ^eing probability that the 

molecule has the rotational energy 


Ej„, 


and jnjm ) is given by 






j~e k(t)dt 


(5.17) 


where 


k(t) = (K-l) <^b[r(0)] T^ m p.(0)] b[r( t) ] (t) ]^> ( 5 . 18) 

and is the correlation function of b( r) T^ m (^'*) , E being the total 
number of molecules. If the spacing between the rotational 
levels of the molecule is small, as is the case for molecules 
having large moment of inertia, it can be assumed that 


J jji_ jit jmi ) j(o) (5.19) 

1 0 

The function J ( 0 ) has been evaluated by Bloom and Oppenheim 
using Constant Acceleration Approximation and an expression 
for this corresponding to dipole-dipole interaction is 
obtained by putting p=2 in eq. (3.49) of Chapter III of this 
thesis, as 

j(0) = j(2) (5.20) 

fi 


where p is the number density of the gas, a is the Lennard- 
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Jones distance parameter, p = 1/kT and \x is the reduced mass 
of the colliding pair and 


l(2) = a-^ s l(2,12) + a2 2 l(2,3) - 2 a^a 2 l ^ 2, 3 y 12) 


(5.21) 


where 1(2,12) etc. are certain integrals as defined in eqs. 
(3-52) and (3.53) in Chapter III. 

Using eq. (4.34) of Rose to evaluate the matrix ele- 
ments occurring in eq. (5.15) and then using eqs. (5.14), one 
gets 


A(JM,J'M') = j(0) [C(J , 1J? M«, M-M')] 2 [ C( J' 1 J,00) ] 2 

( 5 . 22 ) 

The average life time in the state J' is obtained from 

eqs. (5.5) and (5.22) as 

-ijr = ^ 3(0) (5.23) 

T 1 

The spin-lattice relaxation time T^ is then, given by (using 
eqs. (5.3) and (5.2)) 

P T r 2 C 2 

( Ti )“ 1 = 6a *<J(J+1)> (5.24) 


T-, 


or 


8a 2 C e 2 ff 


1 1 Pl a ^ 2lCp|Jr) 

^ <T" > JT~ 


1 


C 2 I ( 2,12) 


+ a 2 2 l(2,3) - 2a 1 a 2 K2,3,12)] 


(5.25) 
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where F = J(J+l), = 2.69 x lO 1 ^ molecules/cc and <> 

represents an average over the rotational, states appreciably 
populated at the temperature of interest. 


2 . Molecules hav ing quadrupole moment 


For this group of molecules (e.g. F^) an expression 
for anisotropic attractive inter molecular potential is 
obtained from eq. (5.9) by putting = SI 2 = 2. The total 
intermolecular anisotropic potential is, then given by 


v A (r,e) = Max) 0(224i m ;L m 2 m) x ) 


a-, 


1 

where b(ax) = 


x' 


12 


a 

x 


2 
3 


and x = r/a and a? = 4mQ 2 /a 5 


(5.26) 


(5.27) 

(5.28) 


where Q is the quadrupole moment of the molecule. 

Proceeding in the same way as for the previous case 
the following relations are obtained 


A( JM$ J'M’ ) 


_ 7 _. 

10m 


0 ( 0 ) 


2£+l 
2 J+l 


[C( J'2J? M' 


M- 


-H* )] 2 [C( J' 2J 5 00)] 2 
(5.29) 



7 

10m 


d(o) [1 - 5 j 2 J+3J-] 


(5.30) 
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T 

and -jj— 
where 

p - ■ 4*34 

i _ jr+ih 
1 5(2J-TTT2J+3T 

1(4) = a-^ 2 I (4,12) + a 2 2 I(4,5) 

V. 3 S ymmetric lo p Molecules 

1 . M olecules hav i ng dipole moments 

Since these molecules (e.g., RH^) have axially symme- 
tric charge distribution, the intermole cular potential 
between two molecules is given by eq. (5.11). The transition 
rates are evaluated in the same way as was done for the 
case V.2.1, the only difference being that the rotational 
states of a symmetric top molecule are described by the 
quantum numbers J, K,M and the factor P JEM is used, in the 
place of P JM » to evaluate the averages. The matrix elements 
are evaluated using normalized symmetric top wave functions 

| JKM> = (|~) * ^ ( opY ) 

and the relation (see eq. (4.30) of Rose.' ) Y pm^ ’ a ^ = ^ 4ri^ ^ 

x u'~ (apQ). After evaluating the matrix elements and carrying 

MO 

out the summations the following expression is obtained for 
the transition rate a( JEM, J'K'M') 


(5.32) 

- 2a 1 a 2 l(4,5 ,12) (5.33) 


8m 2 C 2 ao 

e fi 


7 


lOit 


1 

<F> 


a 


'(2-Jipp.) 2 




1(4) 


(5.31) 
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A(JEM.J'K'M') = (^) j(0) f£±i [C(J'1J ; H',M-M'] 2 

X [ C ( J 1 1 K' s 0 s E) ] 2 6 k , jE (5.34) 
Using eq. (5.6) and then eq. (5.4), t-j ^ and are obtained as 

~J = zfe 3(0) (5.55) 

T 1 

a 1 t 4 N 4 a^( 2 tcP[j.) 2 

^ = 4? U ( 27^ U ■ 

eff 

x [a 1 2 l(2,12) + a 2 2 l(2,3)~2a 1 a 2 l(2,3,12)] ( 5 . 36 ) 

where a = ^ (5.37) 

I being the moment of inertia of the molecule about an axis 
perpendicular to the symmetry axis. 


2 » Molecules ha vi ng quadrupole moment 

For molecules having quadrupole moment and axially 
symmetric charge distribution (e.g. , BF^) the interaction 
potential v.(r,a) is given by eq. (5.26). Evaluating the 
matrix elements etc. in the same way as done for the previous 
case one gets, 

A(Jm,J'K'M') = II 3 - jfj- (|j) 2 3(0) Jj4r 

x [C(J’2Jj M' ,M-M' )] 2 [c( J'2J, K'OK)] 2 6 E , ^ E (5.38) 



■ 84 
•125* 


) 3(0) 


(5.39) 
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L 1 


a 


1 


4tc‘ 


n 2 

°eff 


r 84 1 P 1 z A (2^y.y 

*■ 125t J ^2 I ^ 4 ’ i2 ^ 


(5.40) 


+ a 2 2 l(4 ? 5)-2a 1 a 2 I(4 ? 5,12)3 
where all the quantities are corresponding to the molecules 
of interest. 


V.4 Tq data and Interpre tation 

i'he data as a function of density and temperature 

27 

are available in literature, for hydrogen nuclei in HCu 
and and for the fluorine nuclei in F 2 26 and BF^ 2 and 

are given in table 5.1. The values of the strengths, a ± and 
a 2 ? of the repulsive and attractive terms of the anisotropic 
intermolecular potentials for these molecules were obtained 
by the help of eq.s . (5.25), (5.36), (5.31) and (5.40) respeo- 
tively by the method of least squares, i’he integrals 
I ( 2,12) etc. occurring in these equations were evaluated 
in the manner described in Chapter III and their values are 
given in tables 5.3 and 5.4- The values of other factors i.e., 
c 2 I , a and e/k for these molecules are given in table 
5.2. From the values obtained for a 2 the dipole and quadrupole 
moments were calculated using relations (5.13) and (5.28) 
using appropriate values of the parameter a for the molecule 
under consideration. The values of the electric moments thus 
obtained together with the values reported from other methods 

are given in table 5.5* 
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Tab le 5.1 . Experimental T^/p values for HC1, 

IIH- an d BF V , 


System 


NH. 


5 


bf 3 


HC1 


F, 


apex ux'e 

(°K) 

x l/p 

(msec/amagat ) 

300 

51.3 

348 

40.0 

220 

7.13 

277 

5.43 

300 

4.89 

337 

3. 98 

228 

24.7 

240 

22.2 

263 

19.9 

316 

17.5 

338 

14.8 

222 

1 

O 

rH 

KA 

m 

o 

CO 

H 

290 

12.7 x 10 3 

333 

9.7 x 10~ 3 
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fable 5.2. 

Values of I , 
o ’ 

C eff f £ / ic 8J1<3 - a ^ or 

HC1 etc. 

.System 

I 0 

(1G ^ gm cm 2 ) 

C 2 
ef f 

(MHz) 

s/k S 

(°I0 

a 

U°> 

EC1 

2.71 52 

41x41 55 

360 

3.305 

P 2 

31. 7 26 

26 

157x157^ 

112 

3.653 

BH, 

J 

2.82 25 

217 45 

255 

2.6 

BP 3 

78. 9 25 

42 45 

178 

4.38 


* Values taken from Hirschf elder 


48 
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Table 5.3. Values of integrals 1(2,12), 1(2,3) and 1(2,3 
for a dilute classical gas 


pe 

1(2,12) 

1(2,3) 

1(2,3,12) 

0.73368 

0.02718 

0.20623 

0.05433 

0.85107 

0.02631 

0.21669 

0.05561 

1.06509 

0.02618 

0.23903 

0.05972 

1.13924 

0.02644 

0.24707 

0.06160 

1.36882 

0.02801 

0.27871 

0.06880 

1.50000 

0.02936 

0.29904 

0.07385 

1.57895 

0.03032 

0.31231 

0.07726 


12£ 


Table 5 . 4 . Values of integrals 1(4,12), 1(4,5) and 
1 ( 4 , 5 , 12 ) for a dilute classical gas 


|3£ 


1 ( 4 , 12 ) 1(4,5) 


1(4,5,12) 


0.33630 

0.38620 

0.50450 

0.52819 

0.59333 

0.64260 

0.80909 


0.03125 

0.02805 

0.02345 

0.02288 

0.02156 

0.02088 

0.01947 


0.05024 

0.05009 

0.05082 

0.05111 

0.05205 

0.05295 

0.05678 


0.03363 

0.03193 

0.02967 

0.03572 

0.02943 

0.02879 

0.02909 
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Table 5 

• 5. Values 
for HC1 

of the 

etc. 

electric multi pole moments 


dominant elec- 

Values of the 

moments obtained 

System 

tric multi pole 
moment 

Prom MR 
dat a 

Prom other 
methods 

HC1 

I- 1 - 

(10~ 48 esu 

cm) 

1.25+0.01 

1.18 54 

x 2 

Q 

(10 esu 

cm 2 ) 

0.99+0.01 

0.886 54 

ith. 

1 1 
r 

(10“ 18 esu 

cm) 

0. 29+0.02 

1.47 54 

BP 

3 

Q 

(10~26 esu 

cm 2 ) 

1.32+0.01 

7.99 S 


S A value estimated following the method used by 
P. W. Anderson, Pbys. Rev., 80, 511 (1950). 



150 


It is observed that the value of quadrupole moment Q 
for F2 an< ^ the dipole moment p for HC1 are close to the 
reported values while the values of Q for BP^ and p for 
are quite different. This leads to a conclusion that Rajan's 
model for the correlation time of the intramolecular correla- 
tion function works well in some cases, especially, for sphe- 
rical. top molecules such as CH^ etc. as discussed by them 
and also in linear molecules as found in the present work. 
Perhaps the model does not work well in symmetric top 
molecules. 


45 

It may be noted that Armstrong et al. obtained the 
values of 0^^ from the T^ minimum studies in many symmetric 
top molecules. They found that the values obtained for 
from N.M.R. data are in good agreement with those 
obtained from other methods only for BP^ and CHF^ whereas 
there is a wide descrepancy in the values for CH-^F, PH^ and 
NHy They concluded that the theory for spin -relaxation in 
symmetric top molecules should be reassessed. The present 
analysis also shows a descrepancy only for symmetric top 
molecules. This may be due to the inadequacy of the moael 
proposed by Rajan et al. or due to the inadequacy of the 
spin-relaxation theory for symmetric top molecules itself 
as pointed out by Armstrong et al. 



CHAPTER VI 


CONCLUSIONS 


The work presented in this thesis consists of three 
types of investigations, vizs (1) obtaining information about 
the in ter molecular potentials in the mixtures of H 2 with 
mono-, di- and tri-atomic molecules, namely, Ar, N 2 and N 2 0, 

(2) obtaining the value of spin-rotation coupling constant 
for protons in CH-Cl from the study of T^ minimum, and 

(3) application of a model first proposed by Rajan et al. 
for estimating the correlation time for intramolecular spin- 
rotational interaction. The details of the findings are 
summarized in the following paragraphs : 

(l) Proton spin- lattice relaxation time T-^ in the 

mixtures H 2 ~X, where X is Ar, N 2 or N 2 0, were measured as a 

function of temperature (300-600K for H 2 ~Ar and H 2 ~N 2 and 

200-400K for H 2 -N 2 0) , density and composition with a 30 MHz 

spin-echo spectrometer using phase sensitive detection. A 

Ti—ii / 2—7i pulse sequence was used to produce the echo. The 

T^/ P' versus composition data in' H 2 ~X mixtures were extrapolated 

to 100# of X to find out (T 1 /p-)^*" X , i.e. the proton spin- 

lattice relaxation time due to ortho-H ? and X collisions alone^ 

d as 

by a linear variation of T^/p versus composition as well^by 
using an empirical relation 
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■^ 3 ./ P — Aq + A-^X^. + concentration of X, 

Aq etc. axe constants (6.1) 

to take any non-linearity in 2^/ P versus composition into 

n y 

account. These (l, /p) “ data along -with those in H -Ar and 
H 2 ~N 2 mixtures available in literature at low temperatures 
were analysed using Bloom-Oppenheim theory to obtain informa- 
tion about the intermolecular potentials. The isotropic 
parts of the intermolecular potentials for B^-X Pairs were 
assumed to be 12-6 Lennar d- Jones potential. The anisotropic 
parts were assumed to have repulsive as ire 11 a s attractive 
terms. For H 2 ~h 2 and the radial aid angular dependence 

of the attractive terms were taken to be the same as given by 
the electrostatic multipole interaction between them whereas 
for H 2 - Ar the radial and angle dependences of the attractive 
term were taken appropriate for the electrostatic interaction 
between a quadrupolar linear molecule and a spherically 
symmetric atom. The repulsive terms, for all the three pairs, 
viz. H 2 ~Ar, H 2 ~-N 2 and H 2 -H 2 0 were assumed to vary as r 12 and 
to have the same angle dependence as that of the respective 
attractive terms. The strengths of the repulsive and attrac- 
tive terms of the anisotropic' intermolecular potential were 
obtained by fitting the experimental data. Various forms of 
the anisotropic intermolecular potential were tried to get a 
best representation of the experimental data. It was found 

that the model potentials, written below, for H 0 -Ar, H and 

2 2 ^ 



H 2 -N 2 0 respectively represent the NMR data lest (see eqs. 
(3.81) , (3.84) and (3.85) . 
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[^H 2 .-Ar = 4(e) H 2 -Ar E (f) 12 -(f) 6 ] + 5.59xlO“ 15 erg. 

x [1.45(§) 12 -(”) 6 ]P 2 (cos 9 ^) (6.2) 

with ( e ) jj = 10. 36x10*"'^ org, a = 3.17 R 

f V ^H 2 -N 2 = 4(e) [ (f) 12 -(f) 6 ] + 4.25 x 10" 15 erg 

x [1.64(g) 12 -(f) 6 ][P 2 (cose^) P 2 (cos8> 2 )] 

2 

-27. 73x10 ~ 15 erg (|9 5 ^ a q Y 2q (6 -[) Y 2q (0~p (6.3) 

with ( £ )pi 2 _h 2 = 8*20xlO" 15 erg, a = 3.314 A°. 

C T h,-N,0 = + «-3 x 10- 15 erg 

2 2 d d 

2 

X [2.10(g) 12 -(f) 5 ] ’ l ^a q Y 2q (-h{)Yj q (a 2 ) (6.4) 

with (e) H 0 = 11.54 x1c""'"' erg, a = 3.759 £, 

2 2 

where r is the intermolecular separation, e and a are the 
Lennard-Jones energy and distance parameters for the corres- 
ponding pair, ,0-[ and 6 ! 2 are the angles that the symmetry 
axes of molecules 1 and 2 (colliding pair) respectively make 
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xfith r and and ■ ^ are 3he orientations of the symmetry 

axes -with respect to r. 

Por H 2 ~Ar and pairs the relative anisotropies, 

agAl-e, in the anisotropic attractive r^ term where a 2 is its 
strength, were obtained as 0.137 and 0.131 respectively. 

These are in good agreement with the value (= ; 0.128) obtained 
from polarizability data. This justifies the use of the 
proposed model potentials to explain the NMH data in these 
systems. 

Por H 2 ™P 2 ^ taking the value of the quadrupole moment of 
“-26 

PI 2 as 0.63 x 10 esu cm 2 , the value of quadrupole moment 
of waS obtained from the strength of the attractive 

anisotropic potential as 4.67 x 10 ^ esu cm 2 while the recom- 
mended value is 3.0 x 10~^ esu cm 2 . Hence the proposed 
model potential needs to be improved. 

(2) Proton spin-lattice relaxation times T^ in CH-Cl 
gas were measured at room temperature using the same spectro- 
meter as described in previous paragraphs (a [x/2 burst] -x/2-x 
pulse sequence was used to produce the echo to save the 
experimentation time) and in the density region of 0. 2-6.5 
amagats. The value of 1^ minimum was obtained. By analysing 
the T^/p data the value of the effective spin-rotation coup- 
ling constant was obtained as 0.32 KHz assuming that the 
spin-rotation interaction is the dominent mechanism for proton 
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spin-lattice relaxation in CH^qi gas. The effective cross- 
section for molecular reorientation was obtained as 4.01 if 2 . 
Comparing this cross-section with the kinetic cross-section 
it was found that it takes about 14 collisions to randomize 
the rotational angular momentum, thus the collisions are 
neither very weak nor very strong. 

(3) The existing O^/p data in HC1, f' 2 , and BP^ 

were interpreted using Bloom-Oppenheim theory and a model 
first proposed and used by Raj an et al. which gives an 
estimate for the correlation time of the dominant intramole- 
cular interaction. According to this model the correlation 
time of the spin-rotation interaction, which is the dominant 
mechanism of relaxation in the molecules quoted above, can 
be taken equal to the life time of the molecule in a J 
state. 

The intermole cular potentials for HC1-HC1, Pg-Pg* 

and BP^-BPcj P airs of molecules were assumed to 

consist of isotropic and anisotropic parts, l'he isotropic 

parts were assumed to be a 12-6 Lennard- Jones potential. 

The anisotropic parts were assumed to contain attractive as 

well as repulsive terms. The radial and angular dependences 

of the attractive terms were taken as the electrostatic 

multipole interaction between the molecules whereas the 

-12 

repulsive terms were assumed to have a r radial dependence 
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and the same angle dependence as that of the attractive parts. 

The strengths of the attractive and repulsive parts of the 

anisotropic potentials were calculated from the T^/p data and 

the scaler values of the multi pole moments of the molecules 

were obtained, in turn, from the strength of the corresponding 

attractive terms . The dipole moment p of HC1 and the 

-18 

quadrupole moment Q of F 2 were obtained as 1.25x10 esu cm 

and 0.99xl0 -26 esu cm 2 respectively which are in good agree- 

-~18 

ment with the respective values 1.18x10 ~ esu cm and 
0.886xl0~ 2 ^ esu era 2 obtained by other methods. The values of 
p for and Q for BF^ were obtained as C.29xl0" 18 esu cm 
and 1.32xlO~ 2 ^ esu cm 2 respectively which are far off from 

1 O 

the reported values that are respectively 1.4-7x10 c esu cm 
and 7.99xl0“ 2 ^ esu cm 2 . Thus it may be concluded that the 
model proposed by Raj an et al to find out the correlation 
time of the intramolecular interaction works well in son® 
cases especially for spherical top molecules such as CE^, 

CF 4 and SiF^ as discussed by them and also in linear mole- 
cules as was found in the present work. Perhaps the model 
does not work in symmetric top molecules. 

45 

It may be pointed out that Armstrong et al. obtained 

the value of C 2 from the T, minimum studies in many sym- 
eif u 

metric top molecules. They found that the values obtained 
for C* ff from F.M.R. data are in good agreement with those - 
obtained from other methods only for BF^ and CHF^ whereas 
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there is a wide discrepancy in the values for CH^P, PH^ and 
NH^. They concluded that the theory for spin-relaxation in 
symmetric top molecules should be reassessed. The present 
analysis also shows a discrepancy only for symmetric top 
molecules. This may be due to the inadequacy of the model 
proposed by Raj an et al. or due to the inadequacy of the 
spin-relaxation theory for symmetric top molecules itself 
as pointed out by Armstrong et al. 
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APPENDIX 


CIRCUIT DIAGRAMS 




Pig , A 1 Pulse sequencer 
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Fig. A3 Pulse shaper and amplifier 
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Fig - A6 30 MHz tripler 




Fig . A 7 Power amplifier 








Fig . A 9 Pre - amplifier 





